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Abstract 

This paper is a continuation of [Brlvr] and [Ivr5]. 

Now I start analysis of the Schrodinger operator with the strong magnetic field. 
Keeping in mind the crucial role played by the microlocal canonical form in Chapter 
6 of [Ivrl], devoted to such operators in the smooth case, the lack of the ultimate 
smoothness here seems to be the catastrophic because complete microlocal canonical 
form requires such smoothness. However, combining ideas of section 6 of [Ivrl] and 
[Brlvr] with an idea of the incomplete reduction, with a remainder which is just small 
enough for our needs, I manage under very moderate smoothness conditions to derive 
sharp remainder estimates: namely, under non-degeneracy condition V(V/F) ^ 0, 
assuming a bit more than C 1 and C 2 for d = 3 and d = 2 respectively l \ I recover 
remainder estimates 0(h~ 2 + nh~ l ) and 0(ii~ l h~ l + 1) respectively exactly as in 
the smooth case. Further, for d = 3 without this non-degeneracy condition I recover 
remainder estimate 0(h~ 2 + /i/i" 1 " 1 ^) in C z -case, 1 < I < 2. Here h -C 1 and fi S> 1 
are Planck and coupling constants respectively. 
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Nord. 

1 ^ Actually, the metrics and the magnetic field should belong to C 4 as d = 2 and fi h~ x . 
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Introduction 

0.1 Preliminary 

This paper is a continuation of [Brlvr] and [Ivr5] and a part of the project to derive sharp 
remainder estimates for various classes of irregular operators (i.e. operators with irregular 
coefficients and in irregular domains). 

Here we start analysis of the Schrodinger operator with the strong magnetic field (or, 
equivalently, a large coupling constant /i). 

Keeping in mind the crucial role played by the microlocal canonical form in Chapter 6 
of [Ivrl] , devoted to such operators in the smooth case, the lack of the ultimate smoothness 
here seems to be the catastrophic 2 ) because complete microlocal canonical form requires 
such smoothness. Instead, in this paper I will make a partial reduction with a remainder 
which is just small enough for our needs but not excessively small, and such reduction 
requires much lesser smoothness. 

This reduction serves two purposes. 

First, if magnetic field is not very strong ji < h~a 3 ), I will use it only to prove that 
singularities leave the diagonal and propagate away from it and thus the time interval for 
which we have the standard decomposition of the propagator trace could be extended from 
{|*| <T = efi- 1 } 4 ) to a lar ger one thus improving remainder estimate 0(nh l d ) which is 
due to rescaling as well. 

Namely, for d = 2 I use non-degeneracy condition |V(V/.F)| ^ and extend time to 
{|*| < T\ — efi} thus improving estimate to O^fi" 1 ^ 1 ). For d = 3 no non- degeneracy 
condition at this stage is needed: I consider zone {|^| > C^l" 1 } where £3 is a velocity 
along magnetic field and I use this movement to prove that singularities leave the diagonal, 
increasing time interval to {|t| < e\C,s\} 5 - ) and another zone {|^3| < C/i -1 } which is small 
enough to estimate it contribution to the remainder by 0(h~ 2 ). 

On the other hand, if magnetic field is strong enough, I use this canonical form to 
derive spectral asymptotics as well. For d = 3 and \xh < 1 I do this only in the inner zone 
{IC3I < Cy/JIh}, treating external zone {|^3| > Cy/JIh} in the same way as before. Only at 
this stage qualified estimate of the remainder in canonical form is needed. 

Combining this idea of this "poor man canonical form" with the ideas of section 6 of 
[Ivrl], and [Brlvr], I manage to derive sharp spectral asymptotics under very moderate 
smoothness conditions: namely, under non- degeneracy condition V(V/F) 7^ 0, assuming 

2 ) ln [Ivr2, Ivr3] the irregularities were rare enough in the domain of the really strong field. 

3 ) l skip factors | log h \@ in this preliminary part. 

4 ^The result which is just a rescaling of the standard theory (i.e. the theory with /i = 1). 
^Actually, a bit better but I skip logarithmic factors here. 
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C 1 ' 2 - and C 2,1 -regularity for d = 3 and d = 2 respectively, I recover remainder estimates 
0(h~ 2 + fih" 1 ) and 0(/^ _1 /i _1 + 1) respectively exactly as in the smooth case; for d = 3 
without this non- degeneracy condition I recover remainder estimate 0(h~ 2 + /i/i _1_I +2) in 
C l ' a -case, 1 < I <2 with appropriate a. 

Actually, to recover remainder estimate 0(1) for d = 2 and \i ^> h^ 1 I need to assume 
that the metrics and the magnetic field should belong to C 4 . 

Here and below C l,a denotes a space of functions which |_^J -derivatives are continuous 
with the continuity modulus \x — y\ l ~^ |log \x — y\\ a if / ^ N or I G N and a > 0; for / e N 
and o < one should take / — 1 instead of \l\ ; a will be positive and large enough in most 
of the statements 6 ^. 

In this paper I consider d = 2 and d = 3 only. For d > 4 the sharp remainder estimate 
was not derived even in the smooth case due to lack of the canonical form 7 ''. However, in 
some of the next articles I will to cover multidimensional case as well, using the idea of the 
partial canonical form of this work. However, a higher smoothness will be required. 

I consider Schrodinger operator in the strong magnetic field 

(0.1) A = A + V(x), A =Y, Pj9 jk (x)Pk, 

j,k<d 

Pj = hDj - fiVj{x), n>l 

with real symmetric positive matrix (g^ h ), real- valued Vj, V and small parameter h and 
large parameter \i. 

I assume that the corresponding magnetic field intensity F (x) = ^J2jk v ^ kl (^i^k — 
d Xk Vj) with the standard absolutely skew-symmetric tensor v^ kl 8 -* does not vanish in do- 
main in question. Note that for d = 2 automatically F 1 = F 2 = and thus magnetic field 
is reduced to F 3 . 

For d = 3 we can straighten vector field F = (F l ,F 2 ,F 3 ) and direct it along x% by 
an appropriate choice of variables xi,X2 9 ^ and make V3 = by an appropriate gauge 
transformation A 1— > e~ % ^ h ®( x ) Ae lflh *. I assume that this is has been done and impose 
smoothness conditions to the reduced operator. 

Thus I assume that for this reduced operator 

(0.2) ^3 = 0, d X9 V k = k = 1,2 

6 ) ln contrast to [Brlvr] and [Ivr5] I am not obsessed here to decrease a but I am still pushing hard I 
down. 

7 ) Because of variable multiplicities of the eigenvalues of the matrix magnetic intensity and their "high- 
order resonances" when one of their linear combination with integer coefficients vanishes. 

^With non-zero components equal ±g5 where here and in what follows g = det(g : ' k ) and (gjk) = {g^) 1 - 
9 'To be consistent one should write rather x 1 ,x 2 '. 
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where the second condition is due to the first one and assumption F 1 = F 2 = 0; then one 
can reduce it further to V 2 = and then 

(0.3) F = (0,0,F 3 ), F 3 = g- l *(d xl V 2 -d X2 V 1 ), g = detQ^)" 1 

and the scalar intensity of magnetic field is 

(0-4) F = g£\F 3 \ = (sV 2 - g 12 g 21 )^\(d Xl V 2 - 3^)1, (g jk ) = ^ k )~ l . 

In particular for d = 2 one can assume that g 33 = 1, g 3 ^ — for j — 1, 2; without taking 
absolute value we get pseudo-scalar intensity then. 
Further, I assume that 

(0.5) F > c~ l 

Multiplying operator by F~^ from both sides one can further reduce to the case 

(0.5)* F=l; 

then operator will be J2j,k p i F ~ X 9 ikp k + F~ x (y - h 2 V x ) with 
V\ = F^ J2j k dj(gi k dkF~^ and this perturbation will not affect our results unless d = 2, 
fj, > h~ l (in which case V should be corrected). 

0.2 Smooth case. Survey 

It is known (see Chapter 6, [Ivrl]) that for d = 2 in the smooth case one can reduce operator 
further to microlocal canonical form 

(0.6) /l 2 ^m,n, P ( X 2,^ 1 hD 2 )(x 2 1 +IJ,- 2 h 2 D 2 1 ) m ^- 2n - p h P 

m+n+p>l 

with the principal part 

(0.6)o (V 2 xl + h 2 D\) + V(x 2 , ^ l hD 2 ) 

and for d = 3 in the smooth case one can reduce operator further to microlocal canonical 
form 

(0.7) v 2 

~ l hD 2 )x 

m+n+p+q>l 

x (xj + n- 2 h 2 D 2 1 ) m (hD 3 ) q fj,- 2m ~ 2q ~ 2k - p h p 
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with the principal part 

(0.7) (Ai + h 2 D\) + h 2 D\ + V(x 2 , x 3 , ^ l hD 2 ). 

Decomposing into Hermitian functions of ^h~^x\ one can replace harmonic oscillator 
(h 2 D\ + fi 2 x 2 ) by its eigenvalues (2j + l)fih, j G bZ + . So, basically for d = 2 our operator 
is a kind of mixture of the usual second order operator and a family of /i _1 /i-pdo's while for 
d — 3 it is a kind of mixture of usual second order operator and a family of /i-differential 
1-dimensional Schrodinger operators with respect to 23 which are also yU -1 /i-pdo's with 
respect to x 2 . This explains major differences between d = 2 and d = 3: 

For d = 2 there are sharp Landau levels (2j + l)/j,h (the eigenvalues of the infinite multi- 
plicity as g lk = V = const); without non- degeneracy condition one cannot expect remainder 
estimate better than 0(yu/i _1 ) and with non-degeneracy condition one can improve it up to 
0(/x _1 /i _1 + 1) because only x (/x/i) -1 of Landau levels are not classically forbidden. 

For d = 3 Landau levels are rather bottoms of the branches of the continuous spectrum 
and the non-degeneracy condition is not that crucial: the best possible remainder estimate 
is 0{h~ 2 + jih" 1 ) while the worst possible remainder estimate depends on the smoothness 
of the reduced operator (presuming nothing wrong happened during reduction) and is 
0(h~ 2 + /i/i -1_I +2) in the heuristic accordance with [Ivr5]. 

For d > 4 Landau levels are J2i<k<r(2jk + l)fk^h where > 0, ±ifk are non-zero 
eigenvalues of the matrix magnetic intensity F, 2r = rankF and j G Z +r , and one can think 
operator to be similar to 2-dimensional as rank F = d and to 3-dimensional as rank F < d 
but there is no complete canonical form even if rank F and multiplicities of fk are constant 
because Y2k3kf k can van ish f° r some j G Z r at some points (resonances). Because of this 
there no complete sharp results even in the smooth case (see [Dim] for incomplete and 
[Rail, Rai2, Rai3, Rai4, Rai5] for non-sharp results). However I plan to tackle this case in 
few papers ([Ivr6, Ivr7] etc). 

0.3 Non-smooth case: heuristics as d = 2 

In non-smooth case we do not have complete canonical form but we do not need it. Further, 
if magnetic field is not very strong we do not need even logarithmic uncertainty principle 
e > C \/ yr x h\ logTij to justify quantization of ^(22, £2) ( we do n °t quantize it in this 
case). Consider X2-shift for time T in the propagation. Note that it does not exceed c\i~ x T 
and under non-degeneracy condition |^ 2 V| > e the shift is exactly of magnitude y u~ 1 T 
as T < e/i. This shift is "observable" on the "quantum level" and the propagator trace 
vanishes provided logarithmic uncertainty principle /i -1 T x e > C[i~ l h\ log h\ holds; here 
both Xi- and ^-scales are e. In order to be able to take T as small as T = e/i" 1 we need 
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to pick up e = C\xh\ \ogh\ which in addition to semiclassical remainder estimate Cu~ x h~ x 
brings an approximation error 0(h~ 2 e l \ \ogh\~ u ) for coefficients belonging to C l,cr . This 
error will be small for /i < h l ~ 6 and 5 = 2/(1 + 2) which would take care of /i < h l ~ 6 in the 
smooth case but we do not have such luxury now. 

As n grows larger we need to treat our operator as a family of /i _1 /i-pdo's and then we 
need the logarithmic uncertainty principle e > C \V [i~ l h\ \ogh\ which in addition to semi- 
classical remainder estimate Cfi" 1 ^ 1 brings an approximation error 0(h~ 2 e l \ \ogh\~ a + 

yu/i -1 ^! \ogh\~ a ) where we assume that gi k , F G C l ' a . An actual reduction and calculations 
are the very big part of our analysis. 

For the best results we need to pick e = C min(/i/i| log 7i| , \J u.~ l h\ log h\) and the thresh- 
old is fj, — (h\ log h\)~3. Below it magnetic field is considered to be weak and above it strong 
(and we treat these two cases differently). The second threshold is \i — h~ l ; for \i above it 
we modify our conditions and only finite number of Landau levels matter; we call magnetic 
field extra-strong then. 

0.4 Non-smooth case: heuristics as d = 3 

In contrast to 2-dimensinal case, there is a movement along x 3 with the speed |£ 3 | and for 
time T shift with respect to x 3 is x T • |£ 3 | (provided |£ 3 | > C/i -1 and T < |£ 3 |). This 
shift is "observable" on the "quantum level" and the propagator trace vanishes provided 
the logarithmic uncertainty principle T|£ 3 | x |£ 3 | > Ch\ \ogh\ holds; also we need to assume 
that e x |^ 3 | > Ch\ \ogh\. Here scale with respect to £ 3 is e|£ 3 |. Plugging T = T = e/i -1 
we arrive to inequality 

161 > Q = Cmax(/i _1 , \/ nh\ \ogh\) 

and also we need to take e = Ch\ log/i|/|^ 3 |. Then the contribution of the outer zone 
{l&l ^ q} t° the approximation error will be 0(h~ 2 ) under C 1 ' 2 regularity of the coefficients. 

I 1 2 

Further, in this case one can replace T\ = e|£ 3 | by a bit larger T x = e|^ 3 | log |^ 3 | | and then 
the contribution of this zone to the remainder does not exceed Ch~ 2 f T^d^ which does 
not exceed Ch~ 2 for this increased T\. 

On the other hand, the contribution of the inner zone {|£ 3 | < g} to the remainder does 
not exceed C/j,gh~ 2 which in turn does not exceed Ch~ 2 as \i < (h\ log/i|) - 3. In this case 
we take e = Chloghg -1 in the inner zone and then its contribution to the approximation 
error does not exceed Ch~ 2 as well. In fact we will increase this threshold to h~3. 

For fi > h~s we will consider the outer zone in the same way as above but in the inner 
zone {|£ 3 | < p} we will treat our operator as a family of 1-dimensional Scrodinger operators. 
Then one needs to take e = Ch\ log/i|/|£ 3 | in this zone as well. Now theory breaks in two 
cases. 
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• Without non-degeneracy condition. Here we apply in rather straight forward manner 
arguments of section 2 of [Ivr5]. 

• Under non-degeneracy condition. It is still not very useful to use shift with respect 
to X3, X2, £2 because it lands us with e = C/ih\ log h\ which leads to the total approx- 
imation error 0(h~ 2 ) only under excessive smoothness assumptions. The role of this 
condition is more humble: after |^3| 2 -partition with respect to X2, X3, £2 this condition 
assures that the ellipticity condition fails only for elements with the total relative 
measure x min(|£ 3 | 2 ////i, l) and this leads us in the end to both remainder and ap- 
proximation error not exceeding Ch~ 2 . The logarithmic uncertainty principle requires 
condition |^| > C/i5|log/i|3 which can be weaken to \^\ > Cha. Zone {|^3| < Chs 
is small enough and its total contribution does not exceed Ch~ 2 as well. 

Finally, for fi > logh\~ l the outer zone disappears. 

0.5 Canonical form and the final answer in the non-smooth case 

Our canonical form contains important terms beyond its principal part given by (0.6)o or 
(0.7)o for d = 2, 3 respectively. The next term will be W'(x2, \i~ x hD<2) or W'(x2, x$, n~ 1 hD 2 ) 
respectively with W = W — V where in the case of the Euclidean metrics g^ k and constant 
magnetic field F effective potential W is equal (modulo 0(^~ l \ log to an average of 
V along circumference {{x' 2 , x 'z) '■ \ x 2 ~ x i\ 2 + W2 ~ &I 2 = — H~ 2 V(x2, x$, £2), x' 3 = £3}. 
Then in the smooth case W = — j(j,~ 2 VA±V + 0(/i -4 ) with the Laplacian in the plane 
orthogonal to magnetic lines but nothing better than W < fi~ l is possible in C'-case with 
I < 2. It will affect our answer. 

0.6 Assumptions and results as d = 2 

Now let us describe the results. Let quantify our assumptions first: namely we assume that 
(0.8) c" 1 < < c V£:|£| = l, 

jk 

(0.9) /,F6^(5(0,1)), 

(0.10) 3n G Z + : V + (2ra + l)fihF G C h<T 

(0.11) I v (tt) I + min \V + (2n + l)Ffih\ > c~\ 

(0.12) ip E Cq (B(0, ^)) , 0<^<1 
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with all smoothness conditions understood uniformly, which means that corresponding 
norms do not exceed c. 

Then for uh ^> 1 condition (0.10) is fulfilled with unique n = n and we refer to it as 
(0.10 n ); further in this case we need to check condition (0.11) with n = n only. 

On the other hand, for fih < c condition (0.10) holds or fails for all n < c(fih)~ 1 
simultaneously and for /i/i < 1 non- degeneracy condition (0.11) is equivalent to 

(0.11)' |V(^)| > c" 1 . 

Finally, we assume that (I, cr) y (I, a) w > and K is large enough. 

Theorem 0.1. Let d = 2 and A be a self-adjoint in L 2 (X) operator given by (0.1). Let 
conditions (0.2), (0.5)* ,(0.8)) - (0.12) be fulfilled in B(0, 1) C X C R 2 with 1 < / < 3 and 
with a > 1 as 1 = 1. 

(i) Let (I, a) = (/,cr), /i < Ch~2 and let 

(0.13) V^-c" 1 in £(0,1). 

Then estimate 

(0.14) | J^(x,x,0)-£^(x,0)jij(x)dx\< 

Cfi^h- 1 + Co 2 + Ch- 2 ( f ih) l \\ogh\ 1 -' 7 

holds with standard Weyl expression £^(x,t) = j^h~ 2 (r — V(x)) +y /g where here and 
below g = det(gj k )^ 1 and e(x, y, r) is the Schwartz kernel of the spectral projector of some 
approximate operator A n > . 

(ii) Let (h\ log h\)~3 < /j, < Ch^ 1 and (/,cr) = (2,1). Let condition (0.13) be fulfilled. 
Then for two framing approximations as in footnote n ) estimate 

(0.15) | J (e(x,x,0) -£™ w (x,0) - $™{x))l>{x)dx\ < 

Cfi^h' 1 + C^h^h^l log/i|5-' J 

10 ) Which means that either I < I or I = I, a < a. 

n )One can take two copies A ± of approximation such that (0.14) holds for both and A~ < A < A + in 
the operator sense. We call them framing approximations . This remark holds for theorems 0.1, 0.3, 0.4, 
0.6 in full. 
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holds with magnetic Weyl expression 

(0.16) £ 2 MW (x, r) = J- ]T 6»(r - V - (2n + lfrhF^^Fy/g 



2tt 

n>0 



where 6(t) = r+ Heaviside function and with correction term 

(0-17) £22&0 = 

— J^^(r - W - (2n+ l)/i/iF) - 0(r - W - (2n + l)^hF)^h~ l F+ 

n>0 

wii/i appropriate corrected potential W . 

(^m] Let fi > Ch~ l and (l,cr) y (2,1). Then for two framing approximations as in 
footnote n > estimate 

(0.18) | J (e(*,x,0) -£ 2 MW (*,0) ~ £?Z)i>(x)dx\ < 

C + C fih" 1 {fi^ 1 h)^\ log fi\^ a + C fi 2 (fi^ 1 h)^\ log fi\^-~ a 

holds with £^OTr given again by (0.17) 13 ^. 

(iv) On the other hand, in frames of (Hi) with (0.11) replaced by 

(0.21) min \V + (2n + l)F^h\ > c" 1 , 

estimate 

(0.22) \e(x,x,0)-S^ w (x,0)\<C s fi~ s Wx e 5(0, -) 

12 ) Expression for W is rather complicated and is given in section 4 but for Euclidean metrics it is described 
in subsection 0.5; here I want just notice that the correction term is 0(h~ 2 fi~ l \ log/i|~ <T ). 

13 ) However now W = V + W'h 2 with 

(0.19) W = -—K(4n 2 F + 4n + 3) + F* ^ dj(g jk d k F-*) 

and K a curvature associated with metrics g^F -1 (and this correction is not due to the lack of the 
smoothness); 

(0.20) g^F^ 1 = a 2 8jk ==> k = 2Aloga. 

Furher the sum in (0.18) contains only a finite number of terms one can skip term j-{W — V)h~ 2 there. 
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holds. 

Remark 0.2. (i) Replacement of condition (0.5)* by (0.5) affects only correction term and 
it exceeds remainder only for fih > 1. 

(ii) We get the optimal remainder estimate 0(fi~ 1 h~ 1 + 1) if (I, a) = (2, 1); for jxh > 1 
we need also to assume that (/, a) = (4, 2). 

(hi) We have this optimal remainder estimate in statement (i) iff I > 1, fi < Ch~^~ 1 '^ l+1 ^ \ \ogh\ a/ 

(iv) On the other hand, we can skip £%j^ T in statement (ii) and still have the optimal 
remainder estimate iff I > 1, ji > eh" 1 ^ 1 ^ 1 ^ | log h\~ a ^ l ~ l \ Cases (iii), (iv) overlap for 
(I, a) = (3,1). 

(v) Using rescaling technique we will be able in section 9 to get rid of condition (0.13). 
0.7 Assumptions and results: d = 3 

As I mentioned, three-dimensional case is a kind of mixture of three-dimensional and one- 
dimensional Schrodinger operators without magnetic field. Due to this while we derive 
spectral asymptotics with remainder estimate no better than 0(h~ 2 + fih~ l ) we can always 
get some remainder estimate without non- degeneracy (microhyperbolicity) condition (0.11). 
In particular for not very strong magnetic field sharp remainder estimate holds for sure: 

Theorem 0.3. Let d = 3 and A be a self-adjoint in L 2 (X) operator given by (0.1). Let 
conditions (0.2), (0.5)*, (0.8) - (0.10), (0.12), (0.13) be fulfilled in 5(0,1) C X C R 3 with 
(l,*) = (2,1). ' 

(i) Let ji < h~3 1 \ogh\~3 . Then for two framing approximations as in footnote n > 
estimate 

(0.23) | J (e(x,x,0) - £^(x,0)^tp(x)dx\ < Ch~ 2 

holds with standard Weyl expression S^(x,t) = ^/^(t — V(x)) '^y/g. 

(ii) Let (I, a) = (2, 1). Then above statement remains true for [i < h~s. 

For fi > non-degenerate and degenerate cases are different. First I consider non- 
degenerate case: 

Theorem 0.4. Let d = 3 and A be a self-adjoint in L 2 (X) operator given by (0.1). Let 
conditions (0.2), (0.5)*, (0.8) - (0.12) be fulfilled in B(0, 1)cIcR 3 with (1,2) ^ (I, a) ^ 
(2, 1) and n = in condition (0.10). 
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(i) Let fi < Ch 1 , (l,cr) = (I, a) and let condition (0.13) be fulfilled. Then asymptotics 
(0.24) | / (e(x,x,0) -£ 3 MW (x,0)Wx)dx| < Ch~ 2 + Ch^(fih) l+ ^\ log/i| z+ ^ ff 



holds with magnetic Weyl expression 

(0.25) £ 3 MW (x, r) = -L £ ( r _ y _ ( 2n + 1)^F) >/T 2 F. 



471-2 

n>0 



fnj Lei h 3 < fi < Ch 1 , (I, a) = (2,1) and lei condition (0.13) be fulfilled. Then 
asymptotics 



(0,6) |y(e-(^,0)- £3 ™(M)-&0))^|<^ 
holds holds with magnetic Weyl expression (0.25) and with correction term 

(0-27) f%£ t {x) = 

^ J^((r - IU - (2n + - (r - V - (2n+ l>/iF)*)fi/T 1 F- 



4tt 2 

n>0 



(wi - V2)h 



2U-3 



67T 2 



appropriate corrected potential W 14 - ) . 
(Hi) Let h^ 1 < yU and (I, a) = (2, 1) . Then for two framing approximations as in footnote 
n > estimate 

(0.28) | / (e(x,x,0) -£™ w (x,0))iP(x)dx\ < C^ih~ l 



holds with magnetic Weyl expression given by (0.26). 

Remark 0.5. Using rescaling technique we will be able in section 9 to get rid of condition 
(0.13) in the previous and in the next theorems. 

Now we are going to drop condition (0.11) as d — 3. 



14 ) Expression for W is rather complicated and is given in section 4 but for Euclidean metrics it is described 
in subsection 0.5; here I want just notice that the correction term is 
0{h-i^- l \\ogh\-' y ). 
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Theorem 0.6. Let d = 3 and A be a self-adjoint in L 2 (X) operator given by (0.1). Let 
conditions {0.2), (0.5)*, (0.8) - (0.10), (0.12) be fulfilled in B(0, 1) C X C M 3 with (1,2) ^ 
(I, o~) d (2, 1) and n = in condition (0.10). 

(i) Let h~3 < fi < h" 1 and condition (0.13) be fulfilled. Then for two framing approxi- 
mations as in footnote 1X > estimate 

(0.29) | j (e(x,x,0) - £ MW (x,0) - £^(x)^(x)dx\ < 

holds 15 l 
Further, 

(0.30) \£%£\ <C//H/rt|iog/i|- CT + <V^/r 2 |iog/i|-f 

2_ a- 

an<2 can be skipped as either /i > h l + 2 \ \ogh\ l + 2 or I = 2, a > 0. 
(n,) For u. > estimate 

(0.31) | y (e(x,£,0) - £ MW (x, 0))v^)dx| < Cfih^^h | l g/i|-^ 

holds. 

0.8 Plan of the paper 

1. Preliminary analysis. We make few preliminary remarks. 

2. d = 2 and magnetic field is relatively weak. Using "precanonical form" we prove 
statement (i) of theorem 0.1. 

3. d — 2 and magnetic field is relatively strong. Canonical form. We derive "poor man 
canonical form" for u. > Ch~s\ \ogh\~3: the true (complete) canonical form does not 
exist because we lack the infinite smoothness. However precanonical form exists with 
rather small error but we need to plug corrected potential W instead of V and their 
difference is the source of correction term. 

4. d — 2: Proof of Theorem 0.1. We prove statements (ii),(iii) of theorem 0.1 (i) and 
basically finish two-dimensional case. 

15 ^Note that the second term in the right-hand expression prevails as fi > h^ l /^ l+2 ^ \ log C z + 2 ) . 
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5. d — 3 and magnetic field is relatively weak. We analyze the "outer zone" {|^3| > 
C(fih\ logh\) 2 (if (^3, £3) are "free" ( "non-magnetic" ) variables) and prove that 
the contribution of this zone to the remainder is 0(h~ 2 ); then we prove instantly 
theorem 0.3 for fi < h~z | log h\~s . 

6. d = 3 and magnetic field is relatively strong. Canonical form. We derive "poor man 
canonical form" for fi > Ch~a | log /i | " 3 . 

7. "Standard" spectral asymptotics. We analyze the "outer zone" as < fi < Chr 1 . 

8. Spectral asymptotics. We analyze inner zone {|^3| < C(fih\ log h\)^} and finish the 
proof of theorem 0.1 (for fi < h~s now) and prove theorems 0.4, 0.6. 

9. Spectral asymptotics. Vanishing V or F. We get rid of condition (0.13) and discuss 
consequences of violation of (0.5). 

1 Preliminary analysis 
1.1 Standard results rescaled 

By standard results we mean results for — While rescaling final spectral asymptotics 
by x 1— > fix, h 1— > fih we would get the standard remainder estimate 0(h 1 ~ d ) converted 
into 0((/x/i) 1_d • = 0(fih 1 ~ d ), something which we want to improve, rescaling of the 
intermediate results will be much more fruitful. 

Let us consider some (7, p)-admissible partition in (x, £) and e-mollification with respect 
to x with e depending on (x, such that 16) 

(1.1) C sP - 1 h\\ogh\<e< 1 . 

At first we assume that e is a numerical small parameter. 

Let A be such approximate operator and e(x, y, r) be its spectral projector. Let u(x, y, t) 
be the Schwartz kernel of the corresponding propagator U(t) = e~ lh ltA (we will never 
consider propagator for non-mollified operator). 

Proposition 1.1. Let e be a small numerical parameter and 

(1.2) n < C^h-^loghl- 1 . 

le 'In fact both 7 and especially p could be vectors and then (1.1) should hold for components. 
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(%) Then 

(1.3) | J $ c {r)(d T e{x,x,T) - F t ^ h -i T (xT{t)u{x,x,t))dr^\ < Ch s 

provided 0(r) = 1 for r < —1, 0(r) = for r > —\, \(t) = 1 for \t\ < |, x(t) — f or 
\t\ > 1 17 \ T = Tq = f eo/i _1 wit/i small enough constant > 0, ( > C s /x/i| log /i| 18 ) and 

Xt® = X(t/T), Mr) = Kr/0; 
(ii) Moreover 

(1.5) \F t ^ h -i T (x T {t) -Xf(t))u(x,x,t)\ < Ch s 
with T = Tq and T = f C s h\ log h\, 

(1.6) \F t ^ h -i T Xf(t) (u(x, x, t) - u°(x, x, t)) | < C/i 2 - d 

(1.7) «°(a;, x, t) = j e ih - ltT d T £ MW (x, r) 
and 

(1.8) \F t ^ h -i rX f(t) (u(x, x, t) - u° st (x, x, t)) | < C^ 2 h 3 - d + Ch l ~ d #{h) 
with 

(1.9) U^j^^r) 

where both £ w (x,t) and £ mw (x,t) are constructed for mollified operator A and = 
t'| logt|~ CT . 

17 ^We pick up all auxiliary functions <f>, x etc satisfying assumption of [Brlvr]: 

(1.4) \D a f\ < (CW) (H ~ m)+ Va:\a\<N = e s \\ogh\ 

with large enough exponent m. 
18 -*So C and T also satisfy logarithmic uncertainty principle £ ■ T > C s h\ \ogh\. 



1 PRELIMINARY ANALYSIS 



15 



Proof. Proof of (1.3) is standard because in the scale x \— > fix, t \— > fit, r i— >• fi 1 r, fi \— > 1 
h I— >• fih we have the standard case (with /i = 1). Note that 



1.10) 1 c (r)F^-i T (xTW«(x,x,t))dr = c (/iA)(xr(t)n(x,x,t)) 



t=o 



Inequalities (1.5) and (1.6)-(1.7) also follow from the same rescaling arguments applied 
to standard results of section 5 of [Brlvr]. For (1.8)-(1.9) one should remember how u is 
constructed for \t\ < C s h\ \ogh\ by the successive approximation method; see section 4.2 of 
[Ivrl] and section 5 of [Brlvr]. However, as an unperturbed operator one can also take 

(1.11) A(x, y, hD x ) = J2 9 jk (y)PA + V(y), P j = hD J - fiV^y; x) 

with Vj(y; x) = Vj(y) + (Wj(y), x — y) rather than more standard Vj(y) (where the former 
includes magnetic field and the latter does not) and here we must plug mollified g jfc , Vj, V 
instead of the original (non-mollified) coefficients. 

One can see easily that for such operator the restriction of the spectral kernel to the 
diagonal e(y, y, r) is exactly 8 MW (y, r) and that all other terms of successive approximation 
contribute no more than 0(h~ 2 ). Really, the contribution of the k-th term does not exceed 
Ch~ 2 ~ k T 2h ~ 2 with T = h\\ogh\. Thus contribution of any term with k > 3 is clearly 
0(h~ 2 ) while contribution of the second term is at most 0(h 2 \ \ogh\ 2 ). Analysis of section 
5 of [Brlvr] (based on rescaling technique) shows that one can take T = h in this estimate 
and then the contribution of the second term does not exceed Ch~ 2 as well. It implies 
estimates (1.6)-(1.7). 

To prove (1.8)-(1.9) we need to use the same construction but with the standard un- 
perturbed operator A(y,hD x ) instead; it is given by (1.11) with Pj = hDj — Vj(y). Then 
one can estimate the contribution of the k-th term in the successive approximation by 
C h~ d+k ~ 1 fi k ~ 1 and thus all the terms but the first two contribute less than Cfj 2 h 2 ~ d . The 
contribution of the first term is exactly £ w and the only term we need to consider is the 
second one 

(1.12) (0 + (A - A)G + - G~(A - A)G-)6(x - y) 

where G ± are respectively forward and backward parametrices of evolution operator hD t — 
A. One can see easily that 

A - A = R + R 1 + R 2 , R = (x-y, V y A{y, hD x )) 

where R±, R 2 are second order operators with the coefficients O (/i 2 |x— y| 2 ) and 0($(\x— ?/[)) 
respectively. Then the contributions of the term given by (1.8) with A — A replaced by Rq 
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will be identically due to standard calculations while the similar contributions for Ri, R2 
will not exceed Cfi 2 h 2 ~ d and Ch~ dr d{h) respectively due to the same arguments of section 
5 [Brlvr]. □ 

1.2 Comparing 8 MW and £ MW 

We need to compare £ MW (x, r) and J ijj(x)£ MW (x, r)dx with £ MW (x, r) and J ^(a;)£ MW (a;, r)<ia; 
respectively. The first comparison is easy: 

(1.13) \£ mw (x,t)-£ mw (x,t)\ < 




(1 +nh)h- 3 ti(e) +nh- 2 ti{e)* for d = 3 
(1 +nh)h~ 2 '&(e) + hid = 2 



while the same difference between £ w (x, r) and £ MW (x,r) obviously would not exceed 
Ch- d $(e). 

Proposition 1.2. Under non- degeneracy condition {0.11) 

(1.14) I J tp(x) (^ MW (x, r) - £ MW (x, r))tte| < C(l + ^ih)h~ d d{e) 
for \r\ < e. 

Proof. We need to prove (1.4) only as $(e) < efih; otherwise it will follow from (1.13). 

Consider d = 2 first. Let us introduce 7-admissible partition of unity with 7 = ed(e). 
If \V + (2j + l)fihF — r| > C7 for all j G Z + on some element of the partition, then on 
this element £ MW (x, r) and £ MW (x, r) coincide while on all other elements their difference 
does not exceed C[ih~ x . However, due to (0.11) the total area of the latter elements does 
not exceed 6*7(1 + /i -1 /?. -1 ) which yields the necessary estimate. 

Consider now d = 3. Then 

(1.15) |£ MW (x,r) -£ MW (x,r)| < C(l + ^{xY^h^^e) 

with 7(2) = mm,- \V + (2j + l)fihF\. Due to condition (0.11) function 7(x)~^ is integrable. 

□ 



2 d = 2 and magnetic field is relatively weak 

We start from 2-dimensional case which is more transparent and clean. 
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2.1 Heuristics. I 

The model of such operator is h 2 D\{hD 2 — ^X\) 2 + V(x) and metaplectic map hD\ \— > hDi, 
hD 2 — /iXi i — > — yUXi, Xi i— > Si + fi~ l hD 2 , x 2 ^ x 2 + [i^hDi transforms this operator into 

(2.1) /i 2 /^ + ^2 + y^ Xi + ^-l^ D2j ^ + ^hDt). 

In some sense this operator is modelled by h 2 D\ + /i 2 x 2 + V(fi~ l hD 2 , x 2 ). Then along coor- 
dinates 22,^2 ~ ^~ l hD 2 singularities propagate with the speed and under assumption 
|<%^| > c _1 for time T the shift with respect to x 2 is x /i~ 1 T. To satisfy logarith- 
mic uncertainty principle we should assume that /i~ 1 T • e > Cfi~ 1 h\\og(fi~ 1 h)\ because 
V(n~ l hD 2 ,x 2 ) is yU _1 /i-pseudo-differential operator with respect to x 2 and e is a character- 
istic size in ^-variable. Therefore, Tu = J u(x,x,t)dx is negligible for To < |t| < Ti with 
T == Ce^hl log^" 1 ^)! and Ti = e/x. 

On the other hand, from the standard results rescaled we know Tu for \t\ < T[ = f e/i" 1 . 
In order these intervals to overlap we need e > Cfih\ log/i| and therefore \x < /i" 1 ] log 
So, we pick up the minimal possible value 

(2.2) £ = C/i/i|log/i| 

and then the approximation error will be Ch~ 2 $(fih\ logh\). 

This construction works for pretty large \x\ however, these arguments are optimal only 

for 

(2.3) ^</2 x = /i-3|log/i|-3. 

i 

In the next section we will manage to work with e = C (/i" 1 /^ log(/i~ 1 /i) |) 2 which is less than 

value given by (2.2) if and only if for /i > /2i. We need inequality e > C (/i" 1 /^ log^" 1 ^) |) 2 
even to consider V(fi~ 1 hD 2 , x 2 ) as a legitimate /i -1 /i-pdo. So, in the frames of this section 
V(fi~ l hD 2 , x 2 ) is not even defined properly. The short answer to this obstacle is that we 
should work directly with operator V(x\ + fi~ l hD 2 ,x 2 + /x -1 /YDi) or just not to make any 
reduction but still to be able to reproduce the above arguments. 

2.2 Heuristics. II: Propagation of singularities 

The goal of the next subsection is to prove that at the energy level and below singularities 
propagate with respect to x\, x 2 with the finite speed, and that with respect to "variables" 
Qj x i + Pjk^~ l Pk (with appropriate coefficients (3jk) singularities propagate with 
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the speed not exceeding Co/x x ; further, the same is true for = Qj + Ylki^ 2 @'jkiPkPi 
with appropriate coefficients finally, in the latter case singularities propagate along 
trajectories of the "Liouvillean" field of (V — r)/F 19 \ 

We define j3jk to have Poisson brackets vanish: {Qj, Pk} = mod 0(/i _1 ), or simply 

(2.5) Q 1 =x 1 -F~ 1 P 2 , Q 2 = x 2 + F~ 1 P 1 , {P li P 2 } = -F 
as we would do in the smooth case; then 

(2.6) {A , Qj} = /I" 1 ajkiPkPi; 

hi 

where A = A — V . In the smooth case we would be able to eliminate a.j\.\ modulo ujjA and 
0(/i -2 ) by an appropriate choice of coefficients /3j ki but now a^i have "smoothness" I — 1 
and then f3j ki will have the same smoothness which is not enough to plug them into Poisson 
brackets. However, we can replace coefficients f3'^ ki by their modulo o(l) approximations 
(obtained by mollification with mollification scale which is o(l) but very slowly decaying). 
Then for these new coefficients we have {A°,Q'j} = fi~ 1 ujjA modulo c^/i" 1 ) and 

(2.7) {AQ'j} = fi~ 1 (u J (T + V) + F~ 1 C j V-LU J (r-A)) 

with Cj = (— iy~ 1 d X2 _.. In the smooth case we know from Chapter 6 of [Ivrl] that Uj = 
CjF~ x and so it should be in our case as well: 

(2.8) {A, Q;.} = /i" 1 [Cj (^1) + Uj{r - A)) mod o^" 1 ) 

and this equality will be sufficient to prove that singularities propagate along trajectories 
of (2.4). 

2.3 Propagation of singularities. Rigorous results 

We prove the following statements: 

19 ) Namely: 

(2.4) — = f-i^ 1 C(v- t )/f, with "Liouvillean" field C v d = (— d X2 v, d Xl v), 

denned below 
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Proposition 2.1. 20) . Let d>2, fi < h~ 1+s and let 

(2.9) M> sup 2V/a + £ 

wit/j arbitrarily small constant e > 0. 

Let fa be supported in B(0, 1), fa = 1 in B(0, 2), \ be supported in [—1, 1] (I remind that 
all such auxiliary functions are appropriate in the sense of [Brlvr] ). Let T = Ch\ \ogh\ < 
T<T{ = e . 

Then 

(2.10) \F t ^ h -i rX T(t)(l - fa MT {x - x))faM T {y - x)u{x,y,t)\ < Ch s Vr < e l 
where here and below t\ > is a small enough constant. 

Proposition 2.2. 21 > Let d = 2 and condition (0.5) be fulfilled. 

(i) Let M > sup \FCv/f\ + e with arbitrarily small constant e > 0. Let 

(2.11) T = Ce^hllogh] + C(ph\logh\)* <T<T 1 = e /x. 
Then for /x < /i < 

(2.12) |F^-i t xt(^)(1 - fawn-MQi - x u Q 2 - x 2 )) w x 

x u(x,y,t)(fa }MlJ/ -i T (Q 1 - xt,Q 2 - x 2 )) w \<Ch s Vr, |r| < e t 

where here and below /xq > ; Ci are large enough constants. Here and b w means Weyl 
quantization of symbol b and due to condition (2.11) logarithmic uncertainty principle holds 
and this quantization of the symbols involved is justified, B l means dual operator; 
In particular, for T > C\ this inequality holds for Qi, Q 2 replaced by x±, x 2 ; 

(ii) Let e > be arbitrarily small constant. Let \ be supported in [1 — e\, 1 + ei], where 
€i and fi are as above. 

Let x = y + ji^FCy/p^T with T satisfying (2.11) and T < e±fi. 

Then for [1q < fx < h~ 1+s with large enough constants M (not depending on e) and /io 

(2.13) \F t ^ h -i T x T (t)(l - fa,M^-MQi ~xi,Q 2 -x 2 )) w x 

x u(x, y, t) (fa^riQi ~ Vi, Q2 - V2)) W * I < Ch s Vr, |r| < ei 

In particular, for T > C\ this inequality holds for Qi, Q 2 replaced by x±, x 2 . 

20 ) Finite speed of propagation with respect to x. 

21 ) Propagation along magnetic lines. 
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Proof. Sketch of the proofs Proofs of all three statements are not very different from the 
proof used in Theorem 3.1 of [Brlvr]. Namely, let tp be a function picked up in the proof 
of that theorem as x\ ft was supported in (— oo, 0] and satisfying certain "regularity" con- 
ditions. Then as a main auxiliary symbol we pick up 

(2.14) <p(f(x,t)-e^) with 

f(x,t) = ^(\x-x\ 2 + e 3 )^TM^, 

(2.15) v(f{Qi(x,0,Q2(x,0,t) -e 3 ) with 

f(Q,t) = ^(/i 2 |g - Q\ 2 + e 3 )*T M 1 -, 

(2.16) p(f(Q[(x, 0,Q'2(x, £),*) -e 3 ) with 

/(£', t) = i (/i 2 |g - q - /i-Vt| 2 + 6 3 ) 1 ± e 5 | 

in the proofs of Propositions 2.1, 2.2(i), 2.2(H) respectively where > are small constants. 
We pick signs "=f" analyzing cases ±t > 0. Here £ = FCv if /i _1 T < e 6 in the proof of 

proposition 2.2(h). We extend it to /i _1 T < e by steps. □ 

Proposition 2.3. Let d = 2, I = l,o~ > 1 and conditions (0.5), (0.1 1) be fulfilled. Let \ be 
supported in [— 1, — |] U [~, 1]. Then for 

(2.17) \i*<\t<Chr* 
and T satisfying (2.11) 

(2.18) |iWv(xr(t)r(<Kx)u))| < Ch s Vr : |r| < e. 

Proof. For fi < e'h~\\ log/i|~3 and for ji > e'h~^\ \ogh\~^, /i -1 T > e with small enough 
constant e' this statement follows immediately from proposition 2.2(h). 

However, for fi > e'h~3 | log/i| _ s we need to prove this statement for T £ [Tq,T ]; in 
this case // -1 Tq < e. 

One can assume without any loss of the generality that \Cv/fX\\ > e and P2 = £2; then 
Qi = x\ + a^2 with a = (S^Vi) -1 , where (9 X2 Vi) _1 is disjoint from 0. 
Then we can use function (2.16) with 



(2.19) 



/ = -(/i(g 1 -Qi)±e 5 t). 
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In this case we can take C/i _1 T-scale with respect to x\ and £2 and ChT~ l \ log/i|-scale with 
respect to £2 and £1 provided either u. < e'(h\ \ogh\)~2 or |£? X2 a;| — e' | log /z | — 1 in y with 
sufficiently small e' and then in its e- vicinity (since a > 1). However, one can reach it by 
just rescaling x 2 | — ► ^2 + kx\. □ 

The following corollary follows immediately from proposition 2.2(ii): 

Corollary 2.4. Statement (ii) of proposition 2.2 remains true for all T satisfying (2.11) if 
we redefine x as x = ^r(y) where is Liouvillean flow defined by (2.4)- 

2.4 Proof of Theorem 0.1 (i) 

Now we are ready to prove statement (i) of Theorem 0.1: 
Proof. Due to proposition 2.3 and estimate (1.5) we have 

(2.20) |iWv((xT(*) - XT>(t))T(^{x)u))\ < Ch s 

with T — T x — efi and V = T = Ce~ l h\ \ogh\. 

On the other hand, due to rescaling arguments one can build u(x,x,t) by the standard 

means on interval [-T, T] with T = f efi^ 1 . 

We want these intervals to overlap: T > Tq which is equivalent to inequality e > 
Cfih\ log/i| and we pick up minimal possible e given by (2.1). Then inequality (2.20) holds 
with T = T\ and T' = T. 

Then, since inequality 

(2.21) |F^ A -i T (xT(*)r(^(z)u))| < Ch l - d 

holds for T — T, the same inequality holds for T = 7\ as well. Therefore due to standard 
the Tauberian procedure 

(2.22) |r(e(,.,0)^) -h- 1 J (F^ fc -i T (xT(*)r(«0)))dr| < ^-h l ~ d 

with T — T\ and thus due to (2.20) the same inequality holds with T = T. 
Further, due to proposition 1.1 (ii) we can replace 

h- 1 ! (i^h-ir(xr(*)r«(M)))dr 
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by £ w (x, 0) with an error 0(fi 2 ). Furthermore, the difference between "averaged with 
respect to spatial variables" (multiplied by 4>{x)dx and integrated) functions £ w (x, 0) and 
£ MW (x,0) is 0(/i 2 ) as well. 

Thus semiclassical remainder estimate is 0(yU _1 /i _1 ). Finally, due to proposition 1.2 
approximation errors in "averaged with respect to spatial variables" £ w and £ MW are 
0{h~ 2, d{e)) with e = Cfih\\ogh\ which implies remainder estimate 0(/i _1 /i _1 + fi 2 + 
h~ 2 d(fih\ log /i | ) ) as long as fi < h~ Y \ \ogh\~ 1 . □ 

Remark 2.5. (i) For fi < h~z term /i 2 in the right-hand expression of (0.14) does not 
exceed fi^hr 1 . 

(ii) We can skip 0(fi 2 ) in the remainder estimate if we include in formula the correspond- 
ing term in the standard expression for F t ^h-i T XT{t)u(x,x,t) which is equivalent (modulo 
0(fi 4 h 2 )) to replacement £ w by £ MW . This takes care of fi < h~2 and there is very little 
sense to push it further in this statement. 



3 d = 2 and magnetic field is relatively strong. Canon- 
ical form 

Now we analyze the case 

(3.1) Jii = f h~^\ log/i|~^ < fi < Ji 2 = f h~ l \ log/il" 1 ; 

the second restriction will be removed in subsection 3.4. 

Our arguments will not change much for larger fi as well but some of them become much 
simpler while other become a bit more complicated. We will take 

(3.2) fi' 1 > e > C^ii~ l h\\ogh\ 
which makes it possible to use /i -1 /i-pdos. 

3.1 Canonical form. Main part of operator. I 

Now our goal is to reduce our operator to a some kind of the canonical form. It would 
be easier for constant or at least smooth gi k and F. However, we can overcome technical 
difficulties in a more general well. Instead of normal FIO we will use "poor man 

FIO" . In this section for a sake of simplicity of some arguments we assume that a > but 
one can get rid of this assumption easily. 
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In this subsection we will reduce A to pre-canonical form which is really easy for 
constant metrics g^ k and F and is not needed at all for Euclidean metrics and constant F. 
First of all, we can assume without any loss of generality that 

(3.3) F 3 = F = 1. 

Really, changing orientation, if necessary, we can make F 3 > 0. Then multiplying mollified 
operator by F~2 both from the right and from the left and commuting it with Pj we will 
get operator of the same form with gi k , V replaced by F~ 1 gi k , F _1 V respectively, modulo 
operator h 2 V'(x) where V'(x) is a function, linear with respect to second derivatives of 
gi k and F and thus V does not exceed C(l + ^(e)e~ 2 ); one can see easily that after 
multiplication by h 2 it does not exceed C$(e) which is an approximation error anyway 
(plus Ch 2 which is less than C/i" 1 ^). Surely, such transformation would affect e(x,x,0) 
but we will return to the original operator later. 

In Chapter 6 of [Ivrl] on this step of reduction we applied /i -1 /i-FIO and let us try the 
same now. Namely, let us transform our operator by T(l), where 

(3.4) T (t) = e^ lh ' HL \ L w = ^P^ fc (x)P, 

jk 

with LP* e T l,c . Obviously L w is a Weyl [i~ l h- quantization of 

(3.5) L(x,0 d = i ^ 2 J2L' k (x)p J (x,OPk(x,0, V, =0 - Vj(x). 

In this subsection L, T, t do not denote the same things as everywhere else and J- l,a 
denotes space of functions functions satisfying 

\d:j\<C u {E l -^\\oge\-° + l) 

with the same constants C v as in [Brlvr]. 

Let us consider "Heisenberg evolution" Q(t) = T(—t)QT(t) for /i _1 /i-pdo Q; then 

(3.6) d t Q(t) = i^h" 1 ^, Q{t)\, Q(0) = Q. 
Let us define first symbol q(t) solution to 

(3.7) d t q(t) = {fi~ 2 L, q}, q(0) = q q(t) = q o fa 
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with standard Poison brackets where <pt is a corresponding Hamiltonian flow. Let us consider 
differential equations defining Hamiltonian flow <pt in terms of pj and x k ' 

(3-8) jXj = L jk Pk, j t P> = £ A j k Pk + E /^PfcPm 

A: fc fcm 

with A = f f(x)£J, where A = (A jfc ), £ = (L jfc ) and J = ^ are 2 x 2-matrices, here 

and below & G JF' -1 ^, /(z) = fi^ - 
One can see easily that for |t| < c 



(3.9) \p\ <c =j> |po0 t | <C|p|, 

(3.10) t g 

(3.H) Pj o t = e A ^p + J2^kPjPk on (bl < 

(3- 12) Sj- o<f> t -Xj = ^2 PjkPk 

j 



where p = (puPz), x = (^1,^2)- 
Further, one can see easily that 

(3.13) a o t = 53 + Gfc) = e tA ^ T (g)e tA ^\ 

M t = ^2Pijk;tPiPjPk 

i,j,k 

where a = Ylj kd^PjPk an d (<?)> (<?t) denote corresponding matrices. 

Now, let us consider in this region {\p\ < c/i -1 } third and higher derivatives of apjPk 
with a* G T l ' u ] then 

(3.14) \d^(a PjPk )\ < C v e- V+X \ log^|->- 2 + C 

for 3 < v < c\ \ogh\ where C v are standard constants even if Vj G J rl,cr only. Similarly for 
v > 3 we have 

(3.15) K&{PP&jVk)\ < C u e-» +T - l \ log/x|-V 3 + C; 
so, /i 2 M' G ^fl^-^-V. 
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Now let us pass from symbols to quantizations. I remind that (I, a) >z (2, 0) and also that 
in the standard smooth situations and for a standard quantizations [a w , b w ] = —ih{a, b} w + 
0(h 3 ); then for a standard quantization of /i 2 ao ° <pt equation (3.6) will be fulfilled modulo 
operators with A-bound 22 ' not exceeding C/i 2 (/i _1 /i) 2 x e -1 /! -1 x (e -1 /! -1 + e _2 /i~ 3 ) which 
in turn does not exceed C\T x h. 

Therefore, transformed operator T(—t)A Q T(t) differs from the quantization of the trans- 
formed symbol fi 2 a o <fi t by operator with A-bound not exceeding Cfi~ 1 h. I remind that 
A Q is a Weyl /i _1 /i-quantization of a . 

This is as good as we need. 

Now we can pick up L jfc G T l,a such that e AT (g)e A = ^ g~^j w ^ n 9 = det(<7) -1 

because e A could be any matrix with determinant equal to 1. Therefore we get that modulo 
operator with A-bound not exceeding fi~ 1 h 

(3.17) T(-l)A T(l) = P 2 + P 2 QT 1 ) W ^2 + /U 2 M W , 

/o -i o\ 7i t def , o —12 / dcf , 

(3-18) M = a o ( j) 1 -p 1 - g p 2 , g=go ( p i . 



Z-l,<7 



I remind that according to (3.13) M = £\ . k /3ijkPiPjPk with G .F 
Consider now transformation of V. Note that 

K^Vo^l <C u e l -»\\oghr 

in the same region {\p\ < c/i -1 }. Then one can see easily that modulo operator with A- 
bound not exceeding C(/i _1 /i) 2 xs^ 1 ^" 1 xe l ~ 3 \ log/i|~ CT (which is obviously less than C$(e)) 

(3.19) T(-1)VT(1) = (Vo 0O W . 

Also note that we used only that V\, V 2 G T 1 ^ . Therefore without any loss of generality 
one can assume that 

(3.20) Vi = 0, Pi = hD u V 2 ,d Xl V 2 = g-i e 

where the third assertion follows from the first and the second ones and 

(3.21) P2 = a{x,&)(x 1 -\{x 2 ,&)), a,\eJ*>°. 



22) We say that A bound of B is R if 
(3.16) ||J9u|| < R(\\u\\ + \\A* 

with some exponent m which is fixed in advance. 
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3.2 Canonical form. Main part of operator. II 

In this subsection we will reduce A to canonical form which is really easy for constant gi 
and F . 

Our next transformation is T'(t) = e~ lX " Dl with A w = A w (x2, fi~ 1 hD 2 ); then 

(3.22) T / (-t)DiT'(t) = D 1: T'(-t)A w T'(t) = A w , T'(-t)xiT'(t) = x x + A w 
precisely 

To calculate transformations of other operators we need to introduce the corresponding 
Hamiltonian flow 4>' t : 

(3.23) £i = const, A(x 2 ,^2) = const, x\ = 2i(0) +tA(x 2 ,^2), 
(3-24) jx 2 =ei(%A(x 2 ,6)), ^6 = -6(9, 2 A(x 2 ,6))- 

This flow is less regular than <p t : 
(3.25) \d^' t \ <C7 a Ai- 1 (l + e f - 1 -'l|logAi|- ff ). 

However, using again above arguments we can estimate the difference between operators 
T'(-l)a w T'(l) and {ao^y with a G P> a : its A-bound does not exceed (/i" 1 /*) 2 x fr" 1 /^ 1 x 
e~ 2 fi~ l _ which in turn does not exceed Cfi~ 1 h. The same statement is true for fi~ l /3 w with 
P G T l ~ x ^ . 

Therefore, modulo operator with A -bound not exceeding C^r x h 



(3.26) T , (-1)T(-1)A T(1)T / (1) = A + /i 2 (M') w , 

(3.27) A = P 2 + Pi P 1 = hD 1 , P 2 = -fi Xll 

(3.28) M' = a o $ - a = ^ P' ijk WjVk 

i,j,k 

. . , _ def _ def i def , , / 

with pi = 5, p 2 = - xi, $ = 0i o L . 

Similarly, modulo operator with A Q -bound not exceeding C"d(e) 

(3.29) T'(-1)T(-1)VT(1)T'(1) = (V o $) w . 



3.3 Canonical form. Potential 

Now we have /i _1 /i-pdo A, with the symbol /i 2 a = /i 2 (£ 2 + £?), perturbed by /i 2 (M') w and 
(^/') w with V — (V o $) w . The first perturbing operator would vanish for constant g jfc 
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and F but the second one would still cause trouble because it's symbol depends on 

as well 23 -'. Our goal is to make a perturbation a /i -1 /i-pdo with the symbol depending on 

X2, £2, and r 2 = f x\ + £ 2 only. 

Let us consider the smooth case first as a pilot model; then we can decompose V o $ 
in asymptotic series with respect to To get rid of linear terms we need to make a 

shift in Xi,£i of a magnitude /i -2 ; then the increment of /x 2 ao will absorb these terms. We 
can continue further 24 ) but what we need really is an error in magnetic Weyl expression 
not exceeding C(/i _1 /i + $(e))/i~ 2 . In 3-dimensional case we will need much less than this 
(an error in operator 0(/x -2 ) will be almost sufficiently small) but now we should request 
0(/x _1 /i) error which can be as large as 0(// -4 ) now and to run things in a "smooth" manner 
we would need / = 4 while our objective is to get the best estimate as / = 2 25 \ 

Let us return to the non-smooth case and try to repeat above arguments. To make the 
shift described above we need to use transformation T"(—t) = e~ ltflh 5W with /x -1 /i-pdo S w 
with the symbol belonging to yT 2 T l ' u ] to accommodate M 1 we will need to add /z _1 /i-pdo 
with symbol satisfying improved estimate (3.15). Namely, one can see easily that 

(3.30) \d^M'\ < C v e~ u+1 \ loghl' 1 ^ 3 + C^ +Ul) ~ 

with V\ counting derivatives with respect to (cci,£i) only. 
So in fact 

(3.31) S = S' + S", S' e fT 2 ^, S" e n~ 2 F lfi 
and 5"' satisfies (3.30); some adjustments to be done later. 

Remark 3.1. One can check easily that such operator T"(t) would be /i~ 1 /i-pdo iff 

(3.32) n > m = C7T3|log/i|i 

It is much simpler to analyze transformation by yU -1 /i-pdos than by /i~ 1 /i-FIOs with not 
very regular symbols and one can avoid some hassle if (3.32) holds. This restriction is not 
smoothness related: one just need to take away from fi 2 M' + V a term of magnitude yT x . 

23 )ln this case $ is a linear symplectomorphism. 

24 ^We cannot eliminate terms of even degrees 2m completely, but we can reduce them to /x 2 o; TO (x2, £2X2^1 + 
£ 2 ) m ; see Chapter 6 of [Ivrl]. 

25 ) Actually we need I = 3 because the "weak magnetic field" approach leads to 0(h~ 2 (iJ,h\ log/i|)') error 
and "strong magnetic field" approach in the discussed above smooth version produces 0(/i _2 /i - ') error; 
threshold is as fi = | log h\~^ and to keep error below fi~ 1 h we need / = 3 and a = | and I will proceed 
in this assumption in my next paper [Ivr6] devoted to higher dimensions. 
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One can cover case /i < p\ by results of the previous section; this would require 1 = 1 = 2, 
a = a = 3 to achieve the best possible estimate 0(n~~ l h). 

However I will use more sophisticated arguments to cover fi G [/xi,/x*] here and prove 
the best possible estimate for a = a = 1 instead. 

Let us consider corresponding Hamiltonian flow ipf One can see easily that 

(3.33) \d v x ^t ~I)\< a/x-V" 1 -"! log/i|- CT + C>~ 3 ^ 

where / is an identity map. 

Using arguments of the previous subsection one can prove easily that 
T"(—t)PjT"(t) = fi(pj o ip t ) w modulo operator with the upper A -bound 

(3.34) nifi^h) 2 x (1 + ^"V" 3 ! log h | ^ + xx~ V 2 | \ogh\~ 1 ) x 

(/i-V" 4 | \ogh\-° + /i-V 3 | \ogh\~ 1 ); 

then the same is true for T"(—t)A T"(t) = fi 2 (aooip t ) w . The bad news is that this expression 
is not necessarily less than Cji~ l h even if I = 2, a = 1; the trouble appears as /x is close to 
/i~3 in which case expression (3.34) is close to h. Also, expression (3.34) is not less than 
C'd{e) for I < | and /x close to /i~3 again. 

We could avoid these problems by slightly increasing I and assuming that I > | but 

there is a better way. Note that (3.33) with e replaced by e = f yU _1_<5 with small enough 
<5 > is less than C(/x _1 /i + (with original e in the last expression). Then let us 

assume that 

(3.35) S' G [i~ 2 F l ' a and S" satisfies (3.30) with e replaced by e= f fi~ l ~ s ; here and below 
T* means that in the definition of the class e replaced by e while T* denotes the original 
class. 

The same arguments show that T"(-t)(/i 2 M' + V') w T"(t) = ((/i 2 M' + V) o tp t ) w with the 
same error. So, modulo operator with A-bound not exceeding 
C(/x _1 /i + we need to consider the quantization of 

(3.36) (fi 2 a + n 2 M' + V) ovj) x . 

One can see easily that pj o ip t — pj — t{pj, S' + S"} G /i~ 4 jF z ~ 2,<T and therefore 

/x 2 (a o - ta - {a , S' + 5"}) G iT 2 T l ^ a 
and if V G T l,a and M' satisfied (3.30) with e replaced by e then 

{H 2 W + V) o i/j t - (fi 2 M' + V) would belong to the same class. In this case "the main 
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part" of (3.36) would be /i 2 a + /i 2 {«o, S' + S"} + /i 2 M' + V and we would define S', S" 
from equations 

(3.37) { a ,S'} = -^ 2 (V- W% 

(3.38) {a , = — (M' - fi~ 2 W"). 

Since {a ,g} = (£i<9 xi — x\d^)q these equations would be solvable and solutions would 
satisfy (3.35) if and only if 

(3.39) W'(x 2 ,£ 2 ,r) = ^- J V (r cos t,r sin t;x 2 ,&)dt, 

(3.40) W"(x 2 , r) = ^-// 2 J M'(r cost, r sin t; x 2 , £ 2 ) dt; 

where I remind r 2 = x\ + £ 2 . 

One can calculate easily the difference between (/i 2 ao + fi 2 M' + V) o ipi and /i 2 ao + 
/i 2 {a , 5" + 5"'} + fi 2 M' + V'; it will be 0(fi~ 2 ) which is small enough to be taken care of 
by /i~ 1 /i-pdo transformations in what follows. More precisely 

(3.41) N = (ji 2 a + /i 2 M' + V'^j ofa- 

(/i 2 a + /i 2 {ao, S' + S"} + fi 2 M' + V'^j G /i" 2 ^ 



■1,0- 



Since our symbols belong to T* rather than JF*, we just replace V, M' by their e-mollifications 
V', M' and define S', S", W', W" by (3.37)-(3.40) with V', M' instead of V, M' . 

Then we have unaccounted term K w in the operator before the last transformation with 
K = (/j, 2 M' + V'), M' = M'- M', V' = V- V' and then K G ^(e)^ ' and it also belongs 
to the same class as V + fi 2 M' did. 

To finish this part we need to understand how operator K w is transformed by T"(l). 
One can see easily that K = Kq + K \ + ■ ■ ■ + K m where Kj is the difference between Ej 
and e J+1 -mollifications, Ej = 2 3 e and e m = e. Obviously Kj G r d(Ej)J r ® , ° where T* denotes 
corresponding class with e replaced by Ej. 

Then one can see easily that T"{— 1)KJT"(1) = KJ modulo operator with the norm 
not exceeding [i~ 2 e l ~ 2 \ log/i|~ CT and then T"[— 1)K W T"(1) = K w modulo operator with the 
norm not exceeding fi~ 2 e l ~ 2 \ \ogh\~ a as I < 2 or yU~ 2 | log/i| 1_fT as I = 2 and in both cases 
we got O0d(e)) estimate unless I — 2, \i < Jl\. 

The similar calculations show that for I = 2 and \i G [jtii, jtZ*] 
T" (—1) K W T" (1) = (K') w modulo operator with the norm not exceeding C$(e) with K' = 
K + {K, S' + S"}. 
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3.4 Canonical form. End of reduction 

Thus we arrived to operator A + W w + (L + K') w where W = W + W" , L is defined by 

(3.41) and 

(3.42) K' = K i> K i E c ?Y- 

0<j<m 

Fortunately, both L and K' are small enough to be reduced by fi~ l h-pdo. 

Namely, let us consider transformation by U — (e~ % ^ h ( - s "' +s ) w where S'" and 
s iv 

are symbols of the same class as L and K respectively; in particular, S IV admits 
decomposition of type (3.42). 

One can check easily that e~* M h ( s ""+ 5 ) g JF°'° as long as 
/i^/i^ 1 ! \ogh\~ a < Ce^ 1 which is always the case because \i > Ch~s \ \ogh\~s and (I, a) y 
(1, 2). One can see easily that modulo operator with v4 -bound not exceeding C (fi~ 1 h+ , &(s)) 

(3.43) U~ l (Ao + (W + K' + L) W ^U = A + ({So, S'" + S IV } + W + K + L + R) W 

where one needs to include R G n~ i J-' l ~ 2 ' a only for (/, a) >z (2, 0). 

Now we need just to define S'", S IV , W", W IV from equations similar to (3.37)-(3.40): 



(3.44) {a ,S m } = -(L-W m ), 

(3.45) {a ,S IV } = -(K'-W IV ), 

(3.46) W"'(x 2 ,£ 2 ,r) = — I L(rcost,rsmt;x 2 ,£ 2 )dt, 

2n J 

(3.47) W IV {x 2 , 6, r) = — / tf'(r cos i, r sin x 2 , 6) dt. 

27T ./ 



So, modulo operator with A -bound not exceeding C[yT x h + i?(e) + /i -4 ) we reduced A 
to A + W w where W = W + W" + W" depends on x 2 ,£ 2 ,r only. 

This is almost the end of the story: for (I, a) >z (2,0) and fi G [/2i,/i _ s] one needs to 
make one step more: namely, there is operator R w in the right-hand expression of (3.43) 
and we should make transformation U' = (e - * M h s ) w with S v G u~^T l ~ 2,a defined by 

(3.48) { a ,S"'} = -(R-W v ), 

(3.49) iy y (x 2 , 6, r) = J R{r cos t, r sin t; x 2 , &) dt. 

So, modulo operator with A -bound not exceeding C(u.~ 1 h + $(£)) we reduced A to 
A + Vr w wiere for (/, a) >z (2, 0) one needs to include W v into W as well. 
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3.5 Canonical form. Calculations. I 

Now we will need to calculate W more explicitly. To do this we repeat our construction, 
paying attention to explicit expressions we get rather than their smoothness. Also we will 
consider Weyl symbols rather than operators semselves. Let us start from the beginning of 
subsection 3.3. We started from symbol /x 2 ao + b with b = fi 2 M' + V. Then we applied ipi 
transformation and arrived to 

(3.50) (/i 2 flo + b) o %p l =^ 2 a + 

^ 2 {S, a }+ h 2 {S, {§, a }}+h 2 {S, {S, {S, a }}}+ 
2 o 

b + {S,b} + ^{S, {S, b}} +... 

where S = S' + S" and we skipped terms not exceeding C(/i _1 /i + Plugging 
fi 2 {S, a>o} = — b + Wo where b is a e-mollification of b and 

(3.51) W = M r b = — / b(rcost,rsmt;x 2 ,^2)dt, 

27r Jo 

we get fj, 2 a + W + b\ with 

(3.52) h = Us, -b + Wo} + (6 - 6) + {5, 6} + -Us, {5, -6 + W }} + {5, b}} 

2 b 2 

as a result of subsection 3.3. To calculate W as a result subsection 3.4 (without the last 
step) we need to apply M. r to this expression: 

(3.53) W = W + (M r b) l = (M r b)+ 

M r {S, -k + l -Wo + b} + l -M r {S, {S, ~b + l -Wo + b}}. 

Finally, if / = 2, o > and we want to eliminate 0(/x -4 ) from an error, we need to make 
a last step and define Si from equation (J? {Si, a } = —bi + (W — Wo), calculate 

(3.54) b 2 = \n 2 {S x , {Si, a }} + {Si, Wo} + b x -{W- W ) = 

{Si, -hi + \w+ l -Wo} + bi-(W- Wo) 

and add to W term M. r b 2 . 
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Let W = M r b, b' = | Wo - |6 + b. Then 

(3.55) d^S = ^~ 2 P, with /3 = b-W 0t fi = b-W . 

One can see easily that W G JF Z ' CT and that W = b ( mod 0(#(e) + /i~ 2 )). 
Note that we can take <S = 0asxi=£i = and in this case 

(3.56) \dx*&S\ <Cfi' 2 (fi- 1 + e l - 1 \\ogfi\- a ). 

One can check easily that with the final error O j in all the above calculations 

one can replace standard Poisson brackets {•,•} by "short" ones {-,•}', 
involving only which in polar coordinates (r, (j)) at (x±, £i)-plane are 

{i^fi 1 }' = r_1 ((^/) " (<9rfi0 — (d<f>g) " (d r f))- Then integrating by parts if 6' is differenti- 
ated by <fi we get easily that 

(3.57) M r {S,b'}' = r^drMrib'd^S) = ]^a~ 2 r~ l d r M r {b'^) = 

W 1 = -ji- 2 r' 1 d r M r {{W -bf). 

Further, for (/, a) >z (2, 0) let us consider the last terms in (3.51) and M. r bi 26 \ One can 
see easily that only "short" Poisson brackets should be considered. Then using the same ar- 
guments as above one can prove easily that modulo 
0(jt/ -1 /i + these terms are equal ^yT /L r~ 1 d r M, r {r~' i {b — W-o) 2 ^) and 

^/i~ 2 r~ 1 (9 r .A / l r ((6i — W) 2 ) respectively with b\ redefined as 

(3.52)* bl ^i{ S)b + Wo }, W 1 = M r b 1 . 

3.6 Canonical form. Calculations. II 

Now we want more explicit expressions for W, W\ via g^ k , V. 

First of all, note that for constant metrics g^ h we have b — V — V o where $ : 
T*IR 2 — > T*IR 2 is a linear map. Then Wq = M, r V and since under map $ circle C r = 
{(^2,^2) = const, Xi = r cost,C,i = rsint} translates into ellipse with x-projection S r (y) = 
k9jk( x j ~ Vj)( x k ~ Uk) — r2 } where y = $(x 2 ,£ 2 ) is its center and are elements of 

26 ^One can prove easily that they are not essential for (I, cr) < (2, 0) and we will see that they are x n~ 4 
(l,<r)h (2,0). 



3 D = 2 AND MAGNETIC FIELD IS RELATIVELY STRONG. CANONICAL FORM 33 



inverse matrix (g jfc ) 1 , we get that Wq is M.£ r (^(x 2 , £ 2 )) which is a corresponding average 
of V along £ r (y). 

Similarly, we see that W l = -\^~ 2 \VV'\ 2 = -\^ 2 Y.9 jk {dx V){d Xh V) o $ modulo 
0(/i" 4 ). 

In the general case, the deviation of C r o $ from ellipse S r (with g 3 \, calculated in the 
center) is 0(/i -2 ) and thus W = J^Ag r (^(x 2 ,^ 2 )) modulo 0(/i -2 ). This does not require 
I > 1. It would be nice to get more precise answer in the general well. 

3.7 Canonical form: fi > C 1 h~ 1 \ log h\~ l 

In this case we cannot claim anymore that \pj\ < Cfi -1 in the microlocal sense even if 
fi < ch~ l because f\jip\{x, £), HP2X, £), x) is no more quantizable symbol. However, now 
yU" 1 < £ and f{£~ 1 pi(x,£,),e~ 1 p2(x,£),x) is quantizable symbol and we can claim that 
\pj\ < Ce in the microlocal sense in the zone of r < c + cph due to ellipticity arguments. 
Basically, this makes our construction much simpler; however, we need to remember that 
h 2 is no more non-essential term. 

Construction of subsections 3.1, 3.2 remains the same but we need to take in account 
correction h 2 F~2 Yljk ®j(9 ®kF~* appearing when multiplying A by F~2 from the left and 
right and moving both copies of F~^ to gi k . This correction is the second term in (0.?). 

Further, since the principal part of operator A is of magnitude y,h rather than of magni- 
tude 1 as it was before, the approximation error is 0(/jhe l \ log/i|~°") and it is non-essential 
(i.e. Oifi^h)) for (7, a) = (4,2) only. 

In the very input of subsection 3.2 we can take a Taylor decomposition with respect 
to (xi,£i) (as we did in the smooth case) and then modulo "non-essential symbol" (which 
belongs to (/i" 1 ^ + J 70 ' ) we get 

(3.58) /i 2 (x 2 + £ 2 ) + ^(x 2 ,£ 2 ) + 

3 1 

j=0 j=0 

4 2 

3=0 3=0 

with /3j G J 71 ' 1 , Pj G J 70 ' 1 , ctj G J 71 " 1 **, f3j G J 71 " 2 ** and we need the last term only if I > 2. 

However, we need to remember the terms in decomposition of the main part; we can 
eliminate the third-order term but we will reduce the fourth-order term. More precisely, to 
get rid of the third and the forth terms in (3.58) we use again the construction of subsection 
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3.3 but we look explicitly for 

1 3 

(3.59) S' = J2 6K£i W , S" = ^ 2 J2 4(^2, 6)a&? 1 

j=0 j=0 

with dj, Pj of the same regularity as aij, (3j. 

Now we have expression (3.58) without the third and the fourth terms; the fifth and 
the sixth terms are modified but retain jF' _2 ' <J -regularity provided I = 2, a > 0; for I < 2 
and I = 2, a < we get non-essential term. For /i > c~ l h~ l these terms are non-essential 
as well. So, in these two cases our construction is complete. 

To complete reduction for ji G [h" 1 ] log h\~ l , chr 1 ], I = 2, a > we apply again "poor 
man FIO", this time with 

2 4 

(3.60) S' = &}(*2, G).r|(f S" = v 2 J2 &4tt j 

3=0 j=0 

with &'j,/3j G J rl ~ 2 ' a . However, we cannot eliminate these terms completely, we can only to 
reduce them to 

fe,6)(x?+£i 2 ), «(x 2 ,6)(^ + ei 2 ) 2 

respectively with a,/3 G T l ~ 2 '° '. 

One can calculate easily that the final expression is resembling one from subsection 3.5 
because in our situation M.cf = / + l/-* -2 ^^/ modulo non-essential term and one can 
recalculate this way W\ too. However, construction here reminds one of the smooth theory. 

So, we have proven 

Proposition 3.2. Let fi > h~s | log h\~^ . Then modulo operator with 

A -bound not exceeding C( K n~ l h + $(£)) one can reduce operator A to operator A + W w 

i 

where W = W + W\ + W 2 depends on x 2 , £ 2 and /i^Aq only. 

4 d = 2: Proof of Theorem 0.1 

Now our purpose is to prove statements (ii),(ii) of theorem 0.1. 
4.1 Decomposition 

We again consider more difficult case /i < eh~ l \ \ogh\~ 1 first; the opposite case can be 
considered in the same way with obvious modifications and simplifications. 
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There is very little what is left: removing non-essential operator and compensating 
this by adding =FCo$(e) to the lower/upper approximation (which in turn will not vio- 
late our estimates; see below), we get a 1-dimensional operator- valued /i _1 /i-pdo A + 
W(x2, fi~ 1 hD 2 , fi~ 1 A ). with operator-valued symbol, acting in auxiliary space L 2 (M). 

Taking decomposition into T n (xi) = f (fi~ 1 h)~^v n ({fi~ 1 h)~2Xi) with Hermitian functions 
v n and coefficients which are functions of x 2 we get a family of 1-dimensional /i~ 1 /i-pdos 

(4.1) A n = r 2 n + W(x 2 , ^ l hD 2 , ^-V n ), r n = ((2n + l)fih) * 
as \ihri < Co . 

Even if we cannot apply directly spectral projectors associated with A n we can consider 
propagator e lh lfA and its Schwartz kernel u(x, y, t) and prove easily that for t, \t\ < 7\ = f e/i 

(4.2) u(x, y,t) = T(^2 u n(x 2 , y 2 , i)T n (xi)T n (yi) J T _1 

modulo Schwartz kernel of an operator which becomes negligible multiplied by 
(ijj(x)q(fipi, fip 2 )) where if), q are smooth and supported in B(0, |) and B(0,C ) respec- 
tively; we call operator negligible if its norm does not exceed C(n~ l h) s with large enough 
exponent s. 

Here u is the Schwartz kernel of the propagator e lh tA of the transformed back reduced 
operator A = f t(^A + W(x 2 , fi~ 1 hD 2 , /i~ 1 A ) ^T" 1 27 ) which is 1-dimensional / u~ 1 /i-pdo 

with operator- valued symbol 28 -', T is a full transformation constructed in subsections 3.1, 
3.2, u n are Schwartz kernels of the propagators e lh tAri for A n . Further, by definition 
W = W + Wi because W 2 is not essential in our sense. 

4.2 Remainder estimate 

So, for T = Ti, t < c 

(4.3) F t ^ h - lT (x T (t)T(ij(x)u(x,y,t))^ = 

Ft^h~w(xT(t)^2T(jj(u n (x 2 ,y 2 ,t)T n (x 1 )T n (y 1 )^ J 

^ n>0 ' 

27 ) Apart of just reduction we removed some "non-essential" terms. 

28 ) Actually we have two of them: J 4 ± such that A~ < A < A + and ^4-bounds of operators A — A^ do not 
exceed Ci y ji^ 1 h + ^(e)) (so they are operators described in footnote 10 - ) ). 
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modulo 0(fi~ s ) where as before T is a trace for Schwartz kernel (restriction to the diagonal 
and integration) and ip = T~ 1 tpT is /i -1 /i-pdo rather than a function; its symbol is not 
very regular (see below). 

In this decomposition we can consider only n < C(^ih)^ 1 because for any other n interval 
(— oo, Co) is a classically forbidden zone for operator A n and one can prove easily that their 
contributions do not exceed Cfi~ s (fihn + l)~ s with the arbitrarily large s and thus their 
total contribution is negligible. 

Note that all operators A n are non-degenerate (microhyperbolic) at level in the sense 
that 

(4.4) \An\ + IV^AJ > eo 

if and only if the original operator was: \S7 X {V/F)\ > eo because |V(W — V/F)\ < 
C\x~ x . We assumed that this non- degeneracy condition holds and then in the standard 
way we can prove that if Q x = q(x 2 , /i _1 /i-D X2 ) and Q' y = q'(y2, /x _1 /iZ?j, 2 ) are /i -1 /i-pdos 
with symbols supported in e-vicinity of (x 2 ,£ 2 ) in which condition (4.4) is fulfilled, then 
Ft^h-wXT(t)(QxQ' y Un(x2,y2,t)) is negligible provided 

T < Ti and n~ x T xe> CiT x h\ log/i| or 

(4.5) T = Ce^hl hgh\ < T < Ti = efi; 

I remind that \ is an admissible function supported in [—1, — |]U[|, 1]. Then F t ^ h -i T XT(t)T (ip(x)u(a 
is negligible as well. 

This statement does not look strong enough because T seems to be much larger 
than we need: the direct assault shows that \FijjKu\ does not exceed Ch~ 2 (where K = 
fip(x, £)) W ) is yU _1 /i-pdo cut-off; then 

(4.6) \F t ^ T (x T {t)V{^Ku))\ < Ch~ 2 T 

for x supported in [—1, 1] and equal 1 at [— |, |] and T G [T ,Xi]; here we can take K = I 
because r < c and we can apply ellipticity arguments in zone {\p\ > C'/x -1 }; then (4.6) and 
the standard Tauberian arguments imply spectral remainder estimate 

(4.7) \T(^Q x e)(r) - h~ x £ ^F t ^ T ,(x T (tmQu)))dT'\ < Ch~ 2 ^ 

which is equal to Ch~ 1 fi~ l e~ l \ log/i| and it is close to what we are looking for (but still is 
not that good) only for I = 1. The source of this trouble is that we have e-scale in both x<i 
and £2, while in the differential case treated in [Brlvr] we had scale 1 in £ 2 - 
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However, there is a way to improve our analysis: what we need is an estimate 

(4.8) |iWir(xT(t)r(^))| < C7T 1 + Cfih- 2 -&{e) 

which I am going to prove; then dividing the right-hand expression of (4.8) by T\ we get 
the required remainder estimate 0{yT x Yr l + h~ 2 i9(e)). 
First of all, I want to deal with u n directly. Note that 

i) = ip (x 2l fi~ 1 hD 2 ) + ipi(x 2 , fi~ 1 hD 2 )x l + ip 2 (x 2 , jj,~ 1 hD 2 )/2~ 1 hD l + if>'(x, ^~ l hD) 

with symbol of ip' belonging to fi~ l \ log/i|~°"jF ' ; for / = 2, a > one needs to include 
second order terms 

ipu (%2, yT 1 hD 2 )x\ + ipi 2 (x 2l i~t~ 1 hD 2 )(x l n~ 1 hD l + fi~ 1 hDiXi) + 

iJn(x 2 ,fi- 1 hD 2 )(fi- 1 hD 1 ) 2 

as well. 

Let us consider \F tl _+ h -i T (xt(^)V'' m ) I first; due to analysis above it does not exceed 
Ch~ 2 T x jjT l \ log/i|~°"; one can see easily that it does not exceed the right-hand expression 
of (4.8). 

Plugging terms linear with respect to (x±, fi" 1 hDi) kills trace while plugging quadratic 
terms is equivalent to plugging (n + ^)fi~ 1 hip'(x 2 , fi~ 1 hD 2 ) with if}' = ipn + ip 22 . 

Now we have pure /x _1 /i-pdos, acting on x 2 alone; we consider u n and replace V by V 
acting with respect to x 2 alone. 

Second, we need to improve T slightly, removing | \ogh\ factor and to get better es- 
timates for Fourier transform. Let us consider e-admissible partition in (x 2 ,£, 2 ): 1 = 
^2 u q( u ) and consider element q^) of this partition. Note that for time \t\ < To function 
Q(^(y2,^~ l hD y2 )u is negligible outside Cie-vicinity U( u ) of {x (u)2 , £ (j , )2 ) G suppg^). 

One can easily prove by the standard rescaling technique (see [Ivr5]) the following 
estimate 

(4.9) \T'{ X T(t)q M u n )\ < Ce 2 ^\l^ 

for arbitrarily large s and T e [To, Ti] where 

(4.10) To d =Co£- 1 /i = T /|log/ i |. 
Then 

(4.11) \F t ^ T { X T(t)T\q [v) u n ))\ < Ce 2 ^h- 1 (i^T. 

(4.12) \F t ^ h -i T { X T(t)T'(q {u) u n )) \ < Ce 2 ^- 1 ^ 



4 D = 2: PROOF OF THEOREM 0.1 



38 



for any T G [T ,Ti]. 

Further, due to the standard ellipticity arguments F t ^ h -i T [x-r(t)T' \q( u )U n )) is negligible 
unless 



(4.13) \A n (x 2 ,&)\ < Ce on supp<? ( 



v) 



and for each r, n the set f2 T>n = {(^2,^2) : 1^(^2,^2)1 < Ce} has measure 0(e) due to 
non- degeneracy condition. Therefore we get estimate 

(4.14) \F t ^-i h XT(t)r'(Qu n )\ < Cfih^fo xe = Cfi 
which implies (4.8) and by the standard Tauberian arguments implies 

(4.15) \T(ijQ x e)(r) - h~ l [ (F t ^ h -i Tl ( XT (t)r(il;Qu)))dr'\ < 



1 1 

J. <J11 [ — 

calculate 



with x = 1 on [—hi h], T G [To, T] which is exactly estimate we want. However we need to 



(4.16) h" 1 [ (F^-yfeWrW))^. 

J —00 

Note that we cannot take T = T here, only T = T . 
4.3 Calculations. I 

To calculate (4.16) we rewrite it modulo negligible terms as the sum with respect to n, v of 



(4.17) <S nM d = f 



with q(y) = q^(x 2 ,fJ> 1 hD 2 ). 

One can see easily that on elliptic elements (with (4.13) failed) 

(4.18) S n (y) = ^iihT 1 J J 9(t - A n (x 2 ,C2))'4'(x,C)Q(u)(x2,^2)dxd( 

modulo negligible expression; here 9(r — A n (x2,£2)) = const on suppc^v) and we need to 
calculate the difference between these two expressions on non-elliptic elements (satisfying 
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(4.13)). The total contribution of such terms is not very large: it is less than Ch 1 x/ixex^- 
where the first factor is just present in (4.17), the second one is an estimate (4.14) of 
Fourier transform (F t ^. fi -if l XT(t)T l \Qu n )) , the third is here because this Fourier transform 
is negligible for \r' — r\ > Ce and the last one is the number of n for which (4.13) holds on 
suppg^). Still this expression Ceh~ 2 is larger than what we need. 

We apply method of successive approximations to calculate these terms; as unperturbed 
operator we take 

(4.19) B n = A n (y,^ 1 hD x ), 

fc n = An - B n {y^~ l hD X2 ) = (x 2 - y 2 )Q n (x,y,fi~ 1 hD x ) 

and then 

(4.20) ut(x 2 ,y 2 ,t) =6(±t)u(x 2 ,y 2 ,t) = Tih&(x 2 - y 2 )5(t) = 

=F ihg±6(x 2 - y 2 )5(t) T ihQ±K n g±5{x 2 - y 2 )5(t) 

where Q ± and Q ± are forward/backward parametrices for the Cauchy problem for operators 
hD t — B n and hD t — A n respectively. Similarly 

(4.21) ut(x 2 , y 2 , t)Ql = TthG±5(x 2 - y 2 )Q t y 6(t) = 

T ihQ±5{x 2 - y 2 )QlS(t) T ih^Tl n g^6(x 2 - y 2 )Q t y 6{t) 

for Q y = Q(y 2 , ji^hDyz) where we write pdo acting on y on the right in accordance to ma- 
trix theory and my previous papers, Q l denotes transposed operator (adjoint and complex- 
conjugate). 

So, u n = m° + u\ where u° is a Schwartz kernel of propagator for B n and 

(4.22) v£ = QtTLnut 

We take {x 2 — y 2 ) and carry it forward using commuting relation 

(4.23) [Qt (x 2 - y 2 )\ = Q±[B n , x 2 - y 2 ]Q± 

and in the very front this factor will be killed by T'. Then for our purpose we can replace 

u l n° b y 

(4.24) Qt[Bn,x 2 -y 2 ]QtQnU Q n ± . 
Let us dump it in (4.16). Consider first 

(4-25) F_^- lr L T (t)r'(«(^g^)*)) 
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with either (p? = Xt, T = T Q or ipx = Xt, T > T ; I remind that T = Ch/e. 

In the first case (<pt — Xt, T — To) our estimate is rather straightforward. Note that 
norm of in the strip {\t\ < T} does not exceed /i _1 T, norm of [£> n ,x 2 _ 2/2] does not 
exceed Cfi~ l h. Then factor (x 2 — 2/2) translates into extra factor C pT x h~ x T 2 in the norm 
estimate. For T = T this factor is CpT^he' 2 = C\ \ogh\~ 1 and in the general case it will 



Further, we can assume that norm of TZ n does not exceed Ce l ~ 1 \ logh\~ a + Ce. Really, 
we always can make d X2 W vanish at some point of suppg^) by an appropriate FIO cor- 
responding to some linear symplectomorphism at (X2, £2) -plane. Then in this case (4.25) 
does not exceed Cfie 1 " 1 ] \ogh\~ a ; I remind that for u° it was just Cfi. 

In the second case (y2 T = xt, T > T ) we, using standard rescaling arguments, gain 
extra factor (T Q /T) S with arbitrarily large exponent s which takes care of all (T/T ) k factors 
and summation with respect to partition of unity with respect to t. Now we can conclude 
that (4.25) does not exceed C/ie' -1 ) \ogh\~ a for (p = x, T = T . Then applying arguments 
we used to estimate the total contribution of our "special" elements qi v \ to (4.16), we see 
that the total contribution of these "special" elements qr v ) to (4.16) with u replaced by u 1 
does not exceed Ch~ 2 e l \ \ogh\~ a which is exactly a mollification error estimate. 

Now let us consider contribution of it . One can see easily that it is exactly equal to the 
left-hand expression of 



where T n = T n (xi) and (., .) is an inner product in L 2 {R Xl ). 

The problem is that pg-symbols are not invariant with respect to FIOs and the difference 
between pq- and Weyl symbols is not small because of qi v \. However, there is a simple 
walk-around. Cover our total (x 2 , ^-domain by two, with d&W and d X2 W disjoint from 
respectively. We consider the first domain; the second case will be reduced to the first one 
easily, by means of FIO. 

We can consider partition of unity (depending on n) q'^(x2)q^ l/2 ^2) such that for each 
V\ (and each n) there is no more than one number V2 such that condition (4.13) holds on 



be C\\ogh\- x T 2 lTl. 



(4.26) 





(4.27) 



1 2 



29 )\Vhich means that K = k(x,fj, 1 hD x ). 
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svuppq',, x supp q'/ U2 y Then all the contributions of all other elements do not depend on 
this FIO and summing with respect to v 2 we get just q'/ ui ^(x 2 ). But we can take the required 
linear symplectomorphism (#2, £2) ^ (^2, £2 — r {u 1 )%2) and the corresponding FIO will be 
just multiplication by e~ %llh lr ^) x l/ 2 which does not affect q[ vi -\{x 2 ). On the other hand, for 

ip the difference between different symbols does not exceed CpT x h and therefore without 
increasing our remainder estimate we can rewrite our final answer as 

(4.28) (27r)'V^~ 1 ^ ff 9(T-An(x2,t2))$n(x2,&)dx 2 d&. 



4.4 Calculations. II 

To finish the proof of Theorem 0.1(ii),(iii) we need to bring expression (4.28) with r = to 
more explicit form. 

First of all, we need to use original (xi, x 2 )-coordinates rather than (x 2 , ^-coordinates 
after reduction. The natural expression of the area form in (xi, x 2 )-coordinates is ui = 
y/gdxi A dx2 which can be considered as the area form on E = {pi = p 2 = 0} such that 
ujAdpiAdp 2 = dxAd^. Then after construction of subsections 3.1-3.2 we have a symplectic 
map $ which transforms E into E = {x\ = £1 = 0}; let us denote $|s by <p. Further, $ 
transforms g~?dp 1 A dp 2 at points of E exactly into — dx\ A d£i at points of E . Therefore 
(f)*uj = —dx 2 A d£ 2 . 

Consider next transformations of subsection 3.3. If V = these transformations will 
preserve Eo and be identical on it. In the general case one can check easily that these trans- 
formations transform point (0, 0; x 2 , 6) £ s o info ((i(x 2 , 6)^2(^2, £2); x x + rj 1 (x 2 , £ 2 ),x 2 + 
^2(^2,6)) e E' with |Cj| + |VCj| < C/i~ 2 , \Vj\ + |V?7j| < C/j,~ 4 + Ce l \ loge|~ ,T , and therefore 
after projecting E' onto E : (xi, ^; x 2 , £ 2 ) h_ > (0,0;x 2 ,£ 2 ) we have (1 + k)dx 2 A d£ 2 with 
k = unless (I, a) >z (2, 0) and /i is close to h~s and k = fi~ 4 ki(x) in the latter case. 

So, instead of — dx 2 A d^ 2 in our formula we can use (1 + k) ^Jgdx\dx 2 and instead of 

W(x 2 , £ 2 , p n ) we can use W n = f W(<f)(xi, x 2 ), r n ); further, obviously one can find w(xi,x 2 ) 
of the same regularity as W such that (2n + l)nh + W n < if and only if (2n + l)fih + w < 
and replace W n (x) by w(x). 

One should not be very concerned with the difference between ip and ip n because we are 
interested in calculations when ip = 1 and then ip n = 1 as well. 

Anyway, one can prove easily that 2/^(^2, £2) = 4> 4>~ 1 + Z^ -2 ^ + ^~ Ar >Pn where ?//, 
"0^ are defined by (4.27) with ip replaced by ip', ip" respectively, and with symbols ip', ip" 
belonging to T l ~ x,a , an( j no ^. d e p enc [i n g on ^ h aric [ one can take ip" equal to 

unless (/, a) >z (2, 0) and fi is close to h~ 3 . Then we get the final formula with the correction 
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term 



(4.29) S™, 



—h- 2 j ^2{o(-w - (2n + \)nti)-0(-V - (2n + l)fih))fihil>(x)dx+ 
H^hT 1 / ^29(-w-(2n + l)fih)(iff' n o(j>)dx+ 

J n 

/i" 3 ^" 1 / ^2e{-w - {2n + l)iih){^o<j) + k)dx. 

J n 

The second term is equal to 

/i" 1 ^" 1 J ^Tr! (n^V'^i, x 2 , ix~ 1 hD 1 , £ 2 )) o <^ dx 

where U w is the spectral projector in L 2 (M) for harmonic oscillator [i 2 x\ + h 2 D\ and Tri is 
an operator trace in L 2 {R Xl ). One can apply Weyl formula to calculate this expression. The 
error will not exceed C/i" 1 ^" 1 x (l + s l ~ 2 \ log/i|~ CT ) where the last factor is the remainder 
estimate for Weyl formula for tr(ll w ijj'(xi, hDi)) 30 \ Therefore this error does not exceed 
desired Ch^^h + 0(e)). 

On the other hand, the Weyl answer is 

k- //»(-» - rVt^-Vcos^-Vsina) <* <•*»*.. 

I remind that w — V = O(0(/i -1 )) and then replacing tu by V we make an error not 
exceeding CjjT 2 ~ l \ log h\~ a h~ 2 < C'd{e). We can treat the third term in (4.29) in the same 
way. Therefore 

(4.30) Without deterioration of the remainder estimate the second and third terms in (4.29) 
can be rewritten as Kifi~ 2 h~ 2 and K 2 [i~ A hr 2 respectively with some constants Ki, K 2 . 

I claim that 

(4.31) With an error not exceeding C[i~ l h~ l + Ch~ 2 i9(fih) one can replace in the first term 
in (4.29) the Riemann sum by the corresponding integral. 

30 ) Actually, C is the remainder estimate and Ce l ~ 2 \ \ogh\~ a is an estimate for error in calculations. 
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Really, without any loss of the generality one can assume that | ot,. . V j > eo with either 
j = 1 or j = 2 and the same is true for w. Then the first part of the first term in (4.29) is 
equal to 

(4.32) h~ 2 / ^fihe(r-w-(2n + l)fihYd Xj w){d X] w > )- 1 4j{x)dx+ 

h ~ 2 ^h [ ^0(T-w-(2n+l)/^)(^.™)^ 

J n 

where w' is /i/i-mollification of w. Note that {{d Xj W n )~ l — (d X] W')~ 1 ^ does not exceed 
C(fih) l ~ 1 \ log fih\~ a and therefore passing in the corresponding term from the Riemann sum 
to the corresponding integral brings an error not exceeding Ch~ 2 x fih x (fih) l ~ l \ log /j,h\~ a . 
On the other hand, the first term in (4.32) is equal to 

h~ 2 / ^2(r-w-(2n + l)fih\ {d X] {d X] w')- l ^{x))iihdx 

n 

with (^d Xj (d Xj w')~ 1 ijj(x)^ not exceeding (fih) l ~ 2 \ log/i/i| _<T and passing from the Riemann 
sum to the corresponding integral brings an error not exceeding Ch~ 2 x(fih) 2 x(/j,h) l ~ 2 \ log [ih\~ u 
for (/, a) -< (2, 0); for (/, a) >z (2, 0) we need to repeat this procedure one more time to re- 
cover the same estimate. 

So, (4.31) is proven and with an error not exceeding C{i~ x hr x + Ch~ 2 i9(fih) the first 
term in (4.29) is equal to 

^hr 2 J - V^Jip{x)dx = ^h~ 2 J{V - w)ijj{x)dx 

due to condition (0.13). 

So we get that with indicated error the correction term is 

(4.33) — h~ 2 [(V- w)^{x) dx + h~ 2 (if i/i~ 2 + K 2 ^) 
An J 

with V, Ki,K 2 not depending on /i" 1 and w depending on it. 

Comparing with the results for the weak magnetic field (section 2) when correction 
term is 31 ) we see that expression (4.33) vanishes modulo O^fi^h" 1 + $(e)/i -2 ), with 

31 ^Surely for yu x h~i | log h\~i we can take the same mollification but we also removed operator with 
A -bound not exceeding and one can see easily that it brings error not exceeding Ch~ 2r d(e). 



5 D = 3 AND MAGNETIC FIELD IS RELATIVELY WEAK 



11 



e = (fi 1 h\ log /i | ) 2 and then (4.29) transforms into 
(4.34) £ 2 , corr = 

h~ 2 ( J2{ 9 ( T ~W-(2n+ l)fih)-9(r — V— (2n + l)fih))fihip{x)dx- 

J n 

hr 2 J (V-W)ip(x)dx 

where we replaced notation w by W. 
► Theorem 0.1 (ii) is proven. 

4.5 Calculations: fi > h~ l \ log h\~ l 

In this case calculations are easy. We just should remember the correction h 2 W with W' 
defined by (0.19)-(0.20) to potential. 

5 d = 3 and magnetic field is relatively weak 

In this section we consider the case 

(5.1) fi</2i = C7T3|log/i|-i 

In the smooth case then one needs to reduce operator to canonical form only to ana- 
lyze relatively long (comparing with eo/^ -1 ) propagation of singularities but not to calculate 
asymptotics. The same situation preserves in non-smooth case and to recover sharp asymp- 
totics one needs neither non-degeneracy nor extra smoothness. This is the same threshold 
as for d = 2 despite rather different mechanism. 

5.1 Preliminary remarks 

We consider microlocally mollified Schrodinger operator in the strong magnetic field given 
by (0.1) with symmetric positive matrix (g^), real- valued Vj, V and small parameters 
h, e and large parameter \l. We assume that the corresponding magnetic field F l (x) = 
\ S ^2,jk v '' k \^xj-^-k — d Xk Aj) does not vanish in -B(0, 1); then we can straighten it and direct 
along xi by an appropriate change of variables and make V\ = by an appropriate gauge 
transformation. We assume that this has been done and impose smoothness conditions to 
the reduced operator. 
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Namely we assume that for this reduced operator conditions (0.2), (0.5) hold and 

(5.2) g jk , V, F 3 G .F 1,ff (5(0, 1)), F 1 = F 2 = 0. 
We start from the case (further restrictions to follow) 

(5.3) e > Ch\\ogh\, y, < C- x h\\ogh\- x 
Then 

(5.4) x hr l [A,Xj] = J2(a jkp k + Pk9 jk ){ 

k 

(5.4) 2 h- 1 [A,P 3 ]=h- 1 J2PA9 Jh ,Pi]Pk + h- 1 [V,P 3 ] 

(5.4) 3 ^h-'lA, P t ] = (-l) m F 23 J2(9 5 - l ' k Pk + Pk9 5 - l ' k ) + 

k 

^h' 1 W*> Pi]Pk + ^h' 1 [V, Pi] for z = l,2 

jk 

have symbols of class in any domain {et(x, £) < c/i -2 } of bounded energy where 

(5.5) a(x,0 = a (x,0 + fi~ 2 V, a (x, £) = ^ g 3k Pj(x, £)p k (x, £)> 

Pi(*,0 = &-^)- 

To finish preliminary remarks let us notice that the following statetement holds 

Proposition 5.1. Let / G C^°(IR 3 ) ; / = 1 in 5(0, Co) wt/i /arge enough constant Co. Lei 
T > /i 1 " 5 . Then under condition (5.2) with large enough constant C = C s 

(5.6) \F t ^ h -i T (l - f( Pl (x,t), P2 (x,0,P3(x,0)) W XT(t)u\ <Ch s 

x,ye £?(0,~),r < c. 

One can see easily that / is quantizable under condition (5.3) and the proof is rather 
obvious. 
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5.2 Heuristics 

Let us consider operator (0.1) with Euclidean metrics. Then for time T\ = eo I ^3 1 the shift 
with respect to £ 3 will be less than ||£ 3 | and so £ 3 will be of the same magnitude and 
sign. Then the shift with respect to x 3 will be of magnitude |£ 3 |T because £ 3 is a speed of 
the propagation along x 3 and in order it to be observable we need logarithmic uncertainty 
principle 

(5.7) I&IT-I&I > C/i|log h\ 

because we need to use |£ 3 |-sized scale along £ 3 . This gives us the value T for which we 
can claim that (Tu)(t) is negligible for \t\ G [To,Ti]: To = Ch\ log/i| • |£ 3 |~ 2 . Since we know 

(Tu)(t) for \t\ < T = f e/i _1 due to condition (0.13) and rescaling, we would like to have 
these intervals overlap: T > T or, equivalently, 

(5.8) |£ 3 | > Qi = C^llog h\. 

Then contribution of zone Z{ = f {|^ 3 | > ^i} to the remainder does not exceed J Ch~ 2 f T{ 1 <i£ 3 . 
This integral diverges logarithmically at £ 3 = and is equal to Ch~ 2 \ log Qi\; to improve 
remainder estimate to Ch~ 2 we will need to increase 7\ to e|£ 3 | • [log |£ 3 || CT with a' > 1. 

On the other hand, the contribution of zone Z\ =={|^3| < ^1} to the remainder does 
not exceed Ch~ 2 T~ l Q\ = C[ip\h~ 2 and in order to keep it below Ch~ 2 we need to have 
HQi < C which is exactly condition (5.1). But in this case zone {|£ 3 | < Cji" 1 } is covered by 
these latter arguments we can take Cfi" 1 as a critical value and replace (5.8) by a stronger 
restriction 

(5.9) \^\>Q2 = C^\ 

Condition (5.8) will be restored in sections 6,7 where we analyze the case fi > C(h\ log h\)~3 
and restriction (5.8) is stronger. 

This analysis is basically "smoothness-independent". 

5.3 Variable e 

So far we did not discuss the choice of e; however, since the scale with respect to £ 3 is 
C|£ 3 | in the zone = {|£ 3 | > ^2} the logarithmic uncertainty principle requires e > 
Ch\ \ogh\ ■ l^l" 1 . To reduce approximation error we pick up the smallest allowed e there: 

(5.10) e = Ch\\ogh\ ■ \£ 3 \-\ 
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Then the contribution of Z% to the approximation error will be Ch 3 J $(e)(i6 which is 
0(h~ 2 ) as 1 = 1, a = 2. 

On the other hand, in zone Z 2 = f {|6l < £2} scale with respect to £3 is x p 2 and 
logarithmic uncertainty principle requires e = e 2 == Ch\ log h\ ■ g^ 1 = Cp>{h\ log hi)" 1 (which 
will be also important value in sections 6,7 as it was in sections 3,4). Then the contribution 
of this zone to the approximation error will be Ch~ 3 g~2$(£2) which is 0(h~ 2 ) as I = 1, a > 1. 

However, if we tried to use e = e% everywhere we would get an approximation error 
Ch~ 2q d{ei) which is not 0(h~ 2 ) unless stronger smoothness condition is satisfied. 

So, for d = 3 we will need e to be a function of £3 and may be of x. Really, if we pick 
up the minimal allowed constant e = E2 then approximation error will be 0(h-H{e 2 )) and 
to keep it below h~ 2 we would need to increase I. 

We already met such situation in [Ivr2]. We will pick e to be a temperate function: 

(5.11) /o|V € e|+7|V s e| <ce 
where (7,/?) is (x, £)-scale. In this case we introduce 

(5.12) A = J2^Aai>a 

where tp a are (7, p) -admissible functions, ^* a ^ a is a partition of unity and A a are mollifica- 
tions of A with e = e a calculated at some points of the corresponding elements. 



5.4 Rigorous analysis 

As I already mentioned we consider now the case of the weak magnetic field: p, < fii == |log h\ 
assuming that 

(5.13) 1 = 1, a>l 

and we will not press very hard to reduce a. Let us consider point (x, £) with 

(5.14) g= \Q > Qz^CfT 1 , e = Chg-^loghl 

I remind that the second condition (to e) means logarithmic uncertainty principle; thus 
quantization is possible; we will take e variable but depending on £3 only. 

Let us consider a classical propagation starting from this point. Obviously |^£s| < Co 
due to condition ((5.4) 3 ) which implies that £ 3 keeps it magnitude during time \t\ < e g 
even in the general (not Euclidean) case but we need a bit better. Let us consider 

2 

{/i 2 a(x, 0,6} = P(x, 06 + t* fokPjPk ~ d X3 V 

j,k=i 
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with /3# G JF ,cr . In the smooth case we would be able to eliminate the second term in the 
right-hand expression by replacing £3 by 

2 

(5-15) £ = ^ + J2(3' jkP] Pk 

j,k=i 

with appropriate coefficients (3[ found from the linear algebraic system which we would get 
from equality 

{/i 2 a(x, £),l{x, 0} = ${x, £)£a + fiP"a Q (x, f) - d X3 V mod O^" 1 ) 

(since \pj\ < C/i -1 ). In our case we replace (3'^ k by their /i _<5 -mollifications; then we would 
get the above equations modulo 0(\ log£>|~ CT ). Finally, one can find that [i 2 (3" = d^F and 
therefore 

(5.16) {/i 2 a(x, £),£(x, 0} = F{x, £X(x, + ^ 2 P"a(x, f) - Fd X3 (V/F) 

mod 0(\logg\~ a ). 

Then in an appropriate time direction 32 ) inequality ±£ > \g holds for \t\ G [0, T±] with 

(5.17) T\ = Ti(g) = f eg\ log g\ a . 

This additional logarithmic factor is a small but a crucial progress in comparison with 
T\ = toQ we had before. 

Then using the following function with an appropriate sign q in the standard [Brlvr] 
propagation arguments 

(5.18) <p(±±({l(x,t)-lfat))+os±) 

one can prove the similar statement on microlocal level (due to conditions (5.2), (5 — 4)i_ 3 , 
(5.13)-(5.14) all symbols will be quantizable): 

Proposition 5.2. Let conditions (5.2), (5.3), (5.13) — (5.14) be fulfilled. Let Q be h-pdo 
supported in p-vicinity of (x, £) with |£ 3 | = g > Cfi" 1 . Let Q' be h-pdo with symbol equal to 
1 in g' -vicinity of (x,£) (with g' = Cg\ log g\ a ), intersected with {|£ 3 | > \g}- 

Then (1 — Q')Xl(t)Q is negligible for qt G [0,7]] with 7\ given by (5.17) appropriate 
q — ±1 where here and below U(t) is a propagator associated with A: 

(5.19) U(t) = e~ ih ~ lAt . 

32 ^That means for t of an appropriate sign; namely the sign of —£,3d X3 (V/F). 
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This statement shows that for qt £ [0, Ti] singularity of U(t)Q is confined to ^'-vicinity 
of (x, £) intersected with {|£ 3 | > ^g}. 

Now let us consider classical propagation along x 3 ; then ^x 3 = {fi 2 a, x 3 } = 2g 33 £ 3 + 



/ i Ei=i 2 fl ,3 -'Pi- Replacing x 3 by 

2 

(5.20) 4> = x 3 + ^ ^jPj 

with appropriate coefficients a, found from an appropriate linear algebraic system we would 
able in the smooth case get 

(5.21) {/A, 0} = 2g 3 % mod CK/i" 1 ). 

In the non-smooth case we do the same and no mollification is needed since aj £ JF 1,<T . So, 
the classical shift with respect to <f> for time t with <jt £ [0, Ti] will be at least e gt which is 
"observable" if e g\t\ x g > Ch \ \ogh\. 

On the microlocal level using the same arguments of [Brlvr] with function 

(5-22) ^^(^-^^i^ 

(which again is quantizable due to our conditions) with an appropriate sign <ji = ±1 one 
can easily prove that for qt £ [T ,Ti], 

(5.23) T = Cg~ 2 h\ log h\ 

and ip(x 3 ), supported in e-vicinity of x 3 , ip(x 3 )U(t)Qip(x 3 ) = 0. Then TrXJ (t)Q = 0. Then 
this is true for \t\ £ [Tq,Ti] with no restriction to the sign because of the trace. Therefore 
we have 

Proposition 5.3. Let above assumptions be fulfilled. Then for any T £ [T ,Ti] 

(5.24) iWirUXTxW - XT(*))r(uQ*j) = 0. 

I remind that X: X are admissible functions, supported in [—1, 1] and equal to 0, 1 
respectively on [— |, =]■ 

Then the same statement holds for Q = ip(x)xe(hD 3 ) (we can sum with respect to the 
partition of unity in (x, £1, £2))- From now on operator Q is given by this formula. 



5 D = 3 AND MAGNETIC FIELD IS RELATIVELY WEAK 



50 



On the other hand, from the standard theory rescaled we already know 
Ft-^>h- x T {xt (t)^( u Qy)j with T = eo/i" 1 and in order to combine these two results we need 
to have T > T which is equivalent to (5.8); this condition is wider than our current frames; 
right now we have a/ \ log h\ < which exactly means that /i < /2i = h~s \ \ogh\~a. 

Then from the standard theory rescaled we conclude that 

(5.25) \F t _ hT (xr(*)r(uQ£) ) | < Cgh~ 2 

for any T G [Tq,Ti]. With this estimate the standard Tauberian arguments imply that 

(5.26) K Q ^ \T{eQl){T) - h~ l J* fi^-v (x T (f)r(«Qj)) W| < 

Cgh^T^ 1 = Ch- 2 \logg\~ a . 

Let us take the sum with respect to partition on £3 in zone Z% = {\^3\ > g 2 }- Then the 

right-hand expression of this formula transforms into 

Ch~ 2 J- log g\~ cr dg which does not exceed Ch~ 2 as long as a > 1. 

On the other hand, it follows from the standard theory rescaled that (5.25) holds for 
Q = ip(x)x s (hD 3 ) and T = T = e/i -1 . Therefore for g = g 2 we have TZq < Cg 2 h~ 2 /T = 
Ch~ 2 . 

Combining this estimate with the result of the previous analysis we get the final in- 
equality < Ch~ 2 . 

Further, from the final theory rescaled we know F t ^hr {xT(t)^(uQy)) for T = T and we 
get the final estimate 33 ^ 



(5.27) I J le(x, x, r) - £ w (x, r) J tf>(x) dx\ < Ch~ 2 
where 

(5.28) £ w (x,r) = (27r)- 3 J 6(r - i(z,£)) d£ 

and A is a symbol of the mollified operator and this is h-pdo because we take e depending 
on£ 3 34) . 

33 ^There should be also a term C^hr 1 in the right-hand expression but for fi < pi it is less than Chr 2 . 
I will discuss this expression with more details later. 

34 ) Which means that we take appropriate partition of unity, on each element take its own e, mollify and 
then add. One can check easily that perturbation is small enough to preserve all our arguments. 
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Now let us consider an approximation error in £ . It does not exceed 

(5.29) Ck ^j ^1 log^r^cZ^ = C/i~ 2 | log/i| 1 - CT log/x 

which is 0(h~ 2 ) as a = 2 35 ). 

This finishes the proof of Theorem 0.3(i). 

6 d = 3 and magnetic field is rather strong. Canonical 
form 

From now on we will consider more difficult case 

(6.1) /i> fa = C^/rs | log 

further restrictions will be added later. 

Using arguments of section 5 we will show later (in subsection 7.1.) that the contribution 

of the zone = f {|^ 3 | < Qi}, Q\ = f C{[ih | log h\)^ , to the remainder estimate does not 
exceed Ch~ 2 . 

Unfortunately the trivial estimate in zone Z\ = {\£,3\ < f?i} brings remainder estimate 
0(fiQih~ 2 ) = 0(h~ 2 ^J jj, 3 h\ \ogh\) only. This is not because of the lack of smoothness but 
the deficiency of the method. 

6.1 Canonical form. Preliminary remarks 

To treat this zone Z\ properly one needs to reduce operator to the canonical form as in 
section 6.2 of [Ivrl]. However, as for d = 2, to follow this procedure one must assume that 
s > Cfi~ 1 \logh\ at least and then the approximation error will be too large unless the 
regularity condition holds with I >2. 

On the other hand, these arguments still work, but for much larger fi. One of the 
reasons is that even if we neglect an approximation error, the remainder estimate would 
be at least \xh~ x x 7?-i where is the remainder estimate for 1-dimensional Schrodinger 
operator. Unless we assume that V/F is non degenerate (i.e. V(V/F) ^ 0) we can derive 
nothing better than 1Z\ = 0(h~~~ 1 } according to [Ivr5] and therefore we should not expect 

35 ^Exactly in this place we need to have variable e rather than to take its largest value s(qi) = C/ih\ log h\ 
which would require too restrictive smoothness conditions to get approximation error 0{h~ 2 ). 
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anything better than estimate O^h 1 2 + ; ) for our magnetic 3-dimensional Schr dinger 
operator anyway. For I = 1 this coincides with h~ 2 for fi = h~s exactly. 

However for /i close to we need to modify reduction arguments properly and the 
rest of the section is devoted to this. 

We pay mainly attention to the case 

(6.2) C /r s | log /i| < fi < C^h-^loghl- 1 

and analyze zone Z\. 

First of all we need to decide what does "non-essential" mean. 

(i) Under non-degeneracy condition (0.11) the contribution of zone Z\ to the approxima- 
tion error in £ MW does not exceed Ch~ 3 QiR where R is an upper bound for a perturbation. 

To prove this one needs to notice that for g < g\ A = \/Jih due to this condition the measure of 

the zone y e = f {p 2 < min n > \r — V — (2n+ l)ph\ < 2g} does not exceed Cg 2 / fih and its con- 
tribution to the approximation error in £ MW does not exceed C ' (h~^Qiiih~ 2 g -1 )R x g 2 j \ih. 
Then the total contribution of Z\ to the approximation error does not exceed this expression 
integrated by g~ x dg which gives exactly Ch~ 3 giR. 

On the other hand, under condition (0.11) the remainder estimate will be 0(/i~ 2 ) and 
therefore to keep approximation error below it we need to make calculations with precision 

(6.3) R< R* = Cn-*h?\logh\-*\ 
Therefore in Z\ we can increase e up to 

(6.4) e* = C/i^/iljlog/ill 
as (I, a) = (1,2). 

(ii) In the general case (i.e. without non- degeneracy condition) the contribution of Z\ 
to the approximation error in £ MW does not exceed Ch~ 3 giR + Cph~ 2 R2 while we expect 
the remainder estimate to be 0(h~ 2 + ph~ 1 ~ T +^\\ogh\~~ l ^\ then to keep approximation 
error below it we need to make calculations with precision 

(6.5) R* = Cmin(/i~5/ l l| \ogh\~^ , p~ 2 + h^\ logh\~^) 

where (l,cr) ^ (2, 1). Note that R given by this formula is larger than e l \ \ogh\~ u with e 
defined by (6.1). Therefore in the general case we can increase e up to 

(6.6) e, = C( J R*)T|log/i|T. 
Later further increase will be made. 
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6.2 Canonical form. I 

We assume that g jk G T %1 even if one can weaken conditions. We call term non-essential 
if its v4-bound does not exceed R* given by (6.5); under condition (0.11) we can increase it 
to R*. 

In the zone Z\ we need to have 

(6.7) e = C{fi~ 1 h\logh\)^ 

or larger. Note that for e given by (6.7) mollification error e l \ logh\~ cr does not exceed R*. 

To transform operator A to the canonical form I would like first to transform its principal 
part to the "diagonal form" as it was done for d = 2. This can be achieved by the same 
transformation (3.4) T(t) = f e~ lfl lfl 1<LW with 

L w = \ ^2j k PjL^ k Pk- Then (3.5)-(3.8) are preserved but now we have 3x3 commutator 
matrix J = ^ jj J where J' = K and A = f(x)£J = f^,, jjj with A' = 

/fl2 _/21\ / p tM n\ 

^22 _£2iJ , A" = -^ 31 )- Then e' A = with matrix K(t) = (A')^(e* A - 

7) which is well defined even if A' is degenerate. Therefore we can achieve e AT Ge A to be 
diagonal with the top- left element equal to 1 . However we cannot change g 33 by this method 
unless we use "much larger" operator L (something like u(x)fihD^) which is not a good 
idea at this moment in non-smooth case. 

The rest of the analysis of subsection 3.1 does not change and and therefore we get that 
modulo operator with A-bound not exceeding [i~ x h 

(6.8) T(-1)A T(1) = Pi + P 2 (g') w P2 + WW + A^ w , 

(6.9) M^aoofa-pl-g^pl-g'^pl 

(compare with (3. 17), (3. 18)). 

I remind that according to (3.13) M = YlijkflijkPiPjPk with (3ijk G J 71 ' 1 . Finally, 
modulo non-essential operator T(— l)VT(l) = (V o <pi) w . 

Now, we can assume that g', g" are functions of x only because we can achieve it modi- 
fying M but not changing its type. Further, condition F — 1 means that \fg 1 ] = 1 and we 
can apply the same construction as in subsection 3.1 and therefore we can assume without 
any loss of the generality that 

(6.10) V 1 = V 3 = 0, Pi = hD l =>► V 2 = V 2 (x 2 , x 3 ), 

(6.11) V 2 , d Xl V 2 G J 72,1 
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and 

(6.12) P2 = o;(x,^ 2 )(xi- \(x 2 ,x 3 ,^ 2 )), a, A G J 2,1 . 

Repeating then construction of subsection 3.2, we take the next transformation T'(t) = 
e -iX N D 1 _ ^ s a resu it we w ill get operator 

(6.13) T'(-l)r(-l)A T(l)T'(l) = A + (/i 2 M') w , 

(6.14) A = Pi + P| + P 3 a%P 3 , Pi = hD 1: P 2 = -lixx, P 3 = hD 3 , 

(6.15) M' = a o $ - a = ^ P' ijk pipjpk 

i,j,k 

with pi = f £i, p 2 == — a^i, P3 == £3, $ = f 0i o where 4> t , 4>' t are corresponding Hamiltonian 
flows, «o = 00(^2,^3,^2) £ -T 72,1 (compare with (3.26)-(3.28)). 
Further, modulo operator with Ao-bound not exceeding C'd(e) 

(6.16) T'(-1)T(-1)A T(1)T'(1) = {Vo $) w . 

So far precision was better than we needed: it was e | log/i|~ CT , exactly as in subsections 
3.1, 3.2. We could make ao — 1 but it would be not very useful. 

6.3 Canonical form. II 

Now we will follow subsection 3.3 but with a twist. Namely, we define S f , S", W, W" exactly 
by (3.37)-(3.40) and apply transformation T"{-t) = e^^s™ with phase s = s > + s >t 
Then modulo non-essential term 

(6.17) r\-t)(A + (/i 2 M' + V'Y)T"(t) = ((A + /i 2 M' + F') o^) w 
where ipt is a corresponding Hamiltonian flow. Let us calculate symbol 

(/i 2 a + V 2 M' + V) oip t = ^ao + (/i 2 M' + V) - t/j, 2 {a , S} 

- t{y?M' + V, S} - |{{a , S}, S} + . . . , 

and plug t = 1. Then 

fi 2 {ao,S} = 2^{iid Xl -x 1 d il )S + 2fi 2 a 2 Ud X3 S) + {a 2 , S}$ 
= tf M > + V - Wo + 2 f i 2 a 2 Ud X3 S) + {a 2 , SR 3 2 
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where the last equality is due to (3.37), (3.38) au d W = W + W" . 

Further, let us notice that |£ 3 | < p~ l Qi = C yr x h\ logTij in our zone Z\ where an extra 
factor appears because in the reduction part we consider /i _1 /i-pdo. Then one can see 
easily that modulo non-essential operator (in the general sense) our transformed operator 
equals to the quantization of 

(6.18) /i 2 a + Wo - 2p 2 a 2 &(d X3 S) - ^{/i 2 M' + V - W , S}' 

where {., .}' again mean "short" Poisson brackets (with respect to only). 
Obviously the last term in (6.18) is 0(fi~ 2 ) and it is essential only for 

fi < /ti == eoh~^ | log h\^. Furthermore, d X3 S = O^" 1 "'! log/i|~ CT + /i -3 ) and therefore the 

second term in (6.18) is O ((/i 1 "'] log h\~ a + /i~ 2 ) • log/i|)^) which does not exceed R* 

but is not necessarily smaller than R*. 

Therefore, if condition (0.11) is fulfilled and \i > p,\ then we have finished because only 
the first term in (6.18) is essential. In this case we got operator 

(6.19) Ao + W (x 2 , z 3 , p~ l hD 2 , ^ l hD^ /i _1 A|) w 
with 

(6.20) W (x 2 ,x 3 ,Z 2 ,Z 3 ,p) = M p (fi 2 M' + V). 

Furthermore, under condition (0.11) we will see later that for \i < fii we do not need full 
canonical form at all. 

However, in the general case the second and the third terms in (6.18) are essential and 
we need to continue in the same way as in subsection 3.4. Let us define Si, S'", Wi, W'" 
in the same way as before: 

(6.21) W'" = a 2 d X3 M p S, Wt = ~M p {fi 2 M' + V - W , S}', 

(6.22) 2fi 2 (^d Xl - x 1 d il )S"' = a 2 d X3 (S - M P S), 

(6.23) 2fi 2 (^d Xl - xi%) Si = - 1 -{^ 2 M' + V- Wo, S}' + Wf, 

I remind that a = « (x 2 , x 3 , £ 2 )- Let us apply /x _1 /i-pdo transformation T"'{t) = [e~ %tiih ( S "'^+Si)'\ 
then 

/ 1 \ w 

T"'(-l) (>a + Wo - 2p 2 a 2 Ud X3 S) - -{p 2 M> + V - W , S}') T"{1) = 

(/i 2 a + Wo - 2p 2 alUd X3 S) - l -{p 2 M' + V - W , S}'- 

)W 
■ 
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Plugging (6.21)-(6.22) we see easily that 

(6.24) T w (-l)(/i 2 a + W - 2fjL 2 ^ 3 (d xa S)- 

1 \ w 

~{H 2 M' + V - W , Sy) T"'(l) = 

A + Wx 2 , x 3 , fi~ 1 hD 2 , fi^hD 3 , B) 

where W = W + W'"£ 3 + W x . 

The last step is not really necessary, but we'll do it for the sake of simple canonical 
form. First of all, transformation T v (l), T v (t) = e- ltu)V,D ^ with u = uj(x 2 , x 3 , £2) £ -T 72 ' 1 
transforms D 3 into (d X3 u) w D 3 + ^pt^hDj with (3 = (3(x 2 , x 3 , £2, £3) e ^ r1 ' 1 and there- 
fore one can transform 0^1)3 into D 3 + j3 w fi~ 1 hD 3 . Thus, modulo non-essential opera- 
tor T V (-1)A T V (1) = h 2 D\ + ii 2 x\ + D\. One can see easily that T v (-1)W W T V (1) = 
W w + N w ^- l hD 3 with iV G J^-V. 

Further, modulo non-essential term we can rewrite our reduced operator in the form 
A +W w +N w fj,- 1 hD 3 with a = 1 in expression (6.14) for A and W = W(x 2 , x 3 , £2, V' 1 ^) 
where W = W\ i3=0 G iV = iV(x 2 , x 3 , £2, /^Af), 
iV d i(%lf)| (3=0 6^. 

Finally, applying transformation T y/ (1) with T VI (t) = e~ lt ^ h w and uj' = u'(x 2 , x 3 , £ 2 , ^~ X A^ 
uj' G T ,a we will get (for an appropriate symbol uj') modulo negligible operator Aq + W w . 



6.4 Canonical form. III. fi > eh 1 \logh\ 1 

In this case again an A-bound for n~ l Pj does not exceed Ce and construction is absolutely 
straightforward as in the smooth case; I just note that R* = Ch and 



(6.25) R* = Cmm(h,hw\\ogh\-^) = C 



h for (I, a) di (2,0). 

h\\ogh\~^ for I = 2,0 < a < 1. 



7 "Standard" spectral asymptotics 

In this section I want to prove asymptotics based on F t ^ h -i T {xT^'4 ! Qu) with T = T = e/i -1 
for /j, as large as possible; here Q = q(hD 3 ) is a cut-off operator. First of all, I analyze 

the outer zone Z\ = f \\^ 3 \ > f?i}. Then under condition (0.11) I analyze the inner zone 

2i = {N<^i}. 
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7.1 Outer zone 

First of all, for ji < c _1 /i _1 | \ogh\~ 1 one can apply exactly the same arguments as in section 
5 36) with 

(7.1) q= || 3 | > q x = C(/i/i|log/i|)*, e = ChQ~ l \\ogh\ 

and get an estimate (5.26) for Q with symbol supported in Z\. Then to derive asymptotics 
of T(eQy) (at this stage we include in Q) we need to calculate 

(7.2) h' 1 I' (F t ^ h -i T ,(x T (t)F(uQl))) dr' 
with T = T = eo/i" 1 . Note that 

(7.3) Expression (7.2) with T = T equals modulo 0(h s ) to the same expression with 
T = T = Ch\ log h\. 

Really, let us break Q into few operators Q( v ), with symbols supported in e- vicinity of 
points £(„))■ Without any loss of the generality one can assume that g jk = 5jk at x^) 
37 \ Then if > e (with small enough e > 0) both the propagation speed and scale 

with respect to Xj are disjoint from 0, j = 1, 2, 3 and we can apply arguments of [Brlvr] 
(repeated here many times). 

Now, to calculate (7.2) with T = T we can apply the method of successive approxi- 
mations. However, it is not very convenient to take unperturbed operator in our standard 
way A = A(y,hD x ) because then one can estimate (m + l)-th term of (7.2) with the 
approximation plugged only by 

/ i - 3 |log/i|(/iT 2 /r 1 ) m x /i- 3 |log/i|(/i/i|log/f) m 

which is not good as \i close to h~ l \ \ogh\~ 1 even if we reduce powers of logarithms (which 
is not very difficult). 

Instead we take unperturbed operator 

(7.4) A = J2P^ k (y)h + V(y), P j = hD x .-[iV j (y;x), 

k 

36 >And we need (J, a) = (I, a) = (1,2) only. 

37 )\Ye can achieve it by a linear change of coordinates preserving direction of the magnetic field F. 
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as we did already in 2-dimensional case (1.11). I remind that without any loss of generality 
we assume V 3 = and Vj = Vj(xi,x 2 ) E J- 2 ' 2 . 

With this choice of an unperturbed operator one can estimate (m + l)-th term of (7.2) 
with the approximation plugged by 

h-'T^h- 1 + T 2 h~ l ) m x h-*\ \ogh\ ■ [h\ \ogh\ 2 ) m 

(I remind that \i < Ch\ log h\ ~ l ) . One can reduce the power of logarithms but it is not worth 
of our efforts now; I just notice that modulo 0(h~ 2 ) (7.2) is equal to the same expression 
with u replaced by u° + u 1 where u° and u l are first two terms of approximation: 

(7.5) u° = -ih ^S(x - y)S(t), 

u 1 = -ih sG s 'RJS*8{x - y)5(t) 

where are forward and backward parametrices for (hD t — A) and 1Z = A — A. 

Now we want to get rid of the contribution of u 1 and only extra logarithmic factors in 
its estimate prevents us. Then we need just notice that due to the rescaling method the 
corresponding term in (7.2) with xt replaced by xt with T > T = f Ch~ x does not exceed 
CTh~ A x T 2 h~ l x (-pr) s with an arbitrarily large exponent s. Then the sum with respect 
to t-partition does not exceed Ch~ 2 . On the other hand, this term calculated with xt and 
T = T does not exceed Ch~ 2 as well. Therefore the total contribution of u 1 to the final 
answer does not exceed Chr 2 . 

Finally, let us calculate (7.2) with u = u°. Without any loss of the generality one can 
assume that g^ k (y) = 5jk', we need later to get the answer in the invariant form. Then 

(7.6) u°(x,y,t) = 

h' 1 J j e~ ih ~ X (-'(^ + v (2 ,)) + (s3-^3) diz dT „e y ( Xl , X2] yu y 2] r") 

where e y (xx,X2] Z\,Zz\ t) is the Schwartz kernel of the spectral projector of P 2 + P 2 and u y" 
indicates the same point as in (7.4). Direct calculations (see [Ivrl], Chapter 6 for example) 
show that 

(7.7) e y (x u x 2]yij y 2 ;r") = J2 6 ( T " ~ {^n + l)^hF{y))F{y)^h~ l 

n>0 
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and 

(7.8) T y (u°Ql) = f ^ e ^- 1 *fe 2 + (2n + D^(,) + v(,)) g(e3) ^ 3 x F{y) ^ h -^ 

J n>0 

(7.9) iWv(*r(*)r(uQy) = 

/ f(T/i- 1 (r / -(2n + l)F(y)^-y( 2/ )-e 3 2 ))x 

n>0 

g(&) d& x F(y)^/i- 2 

where x is a Fourier transform of x and modulo 0{h~ 2 ) (7.2) is finally equal to J £q W (x, t)^(x) 
with 

(7.10) ^ w =£/K T - ( 2 " + - ^ - ^teO d& x F/i^ 2 . 

n>0 

One can check easily that e\ log/i|~ 2 < £> 2 for g > Qi and therefore the approximation error 
of £% w does not exceed C J h~ 3 e\ \ogh\~ 2 dg < Chr 2 . So we get 

Proposition 7.1. Let (I, a) = (1,2) and /2i < /i < c^/i" 1 ) log/i| _1 . Lei Q = q(hD 3 ) with 

Q = 1 - Xsi (&) 
T/ien 

(7.11) 1 1 (r.(eQj) -^ w (x,r))^(x)rfx| < C^ 2 
u»& £^ w (x,t) denned by (7.10). 

7.2 Non-degenerate case 

Now I want to prove that under condition (0.11) TiuipQ^) is negligible for T < \t\ < e as 

(7.12) e = C{ih\\ogh\ in Z x 
and Q = q(hD 3 ) with symbol supported in Z\. Namely: 

Proposition 7.2. Let condition (0.11) be fulfilled, Q = q(hD 3 ) with appropriate symbol q 
supported in Z\ . Then 

(7.13) \F t ^ h - lT ( XT r^u)\ < Ch s Vr : |r| < e 
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as T e [To,Ti] where x (in admissible function supported in [— 1, — ~] U [s, 1] and 
(7.14) T = Ce^hl \ogh\, T x = e. 

In particular (7.13) holds as T e [To, e] ratt To = e^^T 1 /or £ given by (7.12). 
Corollary 7.3. Under conditions (0.11), (0.13), (7.12) 



(7.15) | J ( e(x, x, t) - J F t ^ h -i T ,(xfu(x,x,t)) )ip(x)dx\ < Ch~ 

(7.16) | I (e(x, x, r) - £ MW (x, rj)ij>(x)dx\ < ChT 



Remark 7.4. (i) Note that under condition (0.11) and e given by (7.12) the contribution 
of zone Z\ to the approximation error does not exceed 

C(/^|log/i|) /+ 2|log/i|- CT /i- 3 which is in turn does not exceed Ch 2 as lone; as 



(7.17) fi < fjt M = Ch~^\\ogh\- 1+ —i 

while remainder estimate does not exceed Ch~ 2 anyway. 

(ii) Note that /i(i,2) = h~s \ log /i| s and fMn a ) ~ h~~ 1+s in the smooth case when I is large. 

Now statement (i) of theorem 0.4 follows from propositions 7.1, 7.2, corollary 7.3 and 
remark 7.4 

Proof, of Proposition 7.2. Making ex-partition we can assume that on it |VV — i\ < e at 
each point with some fixed £, \£\ > 4e. 

Let us consider first the case when \£±\ > |, £j_ — (t\,t 2 ) and examine propagation of 
singularities with respect to (#2, £2) f° r the transformed operator. 

To time T the shift with respect to x 2 will be x e\i 2 \{i~ x and the logarithmic uncertainty 
principle means that e^i^fi^Txe > C/i _1 /i| log/i| which for \£ 2 \ > e is equivalent to T > T 
with T given by (7.4). 

Similarly, to time T the shift with respect to £2 will be x e|£i|/i _1 and the logarithmic 
uncertainty principle means that e 2 \£i\fi~ 1 T x e > CpT l h\ log h\ which for j^l > e is 
equivalent to T > T again. 

We can justify this analysis by our standard propagation arguments with functions 

(7-18) ^——ie-), ^^ r ±e-) 

respectively. 
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On the other hand, let us consider the case \£ 3 \ > | and examine propagation of sin- 
gularities with respect to £1. To time T the shift will be x e\l 3 \ and and the logarithmic 
uncertainly principle means that e^l/i -1 x e > C7i|log/i| which is equivalent to T > T 
again. 

We can justify this analysis by standard means with a function 

(7.19) <pfkzM± e l.\ 

In both cases upper bound for T is a small constant because for this time (x 2 , x 3 , £ 2 , £3) 
keeps in ei-vicinity of the point of origin 38 ) which can be justified by standard means with 
functions 

(7.20) <p[—±en — ), <^(-±e/i — - — ), 

ft 6 ~V3\ ft x 3 -y 3 , 

^yf ~f~'' ^ t >' 

□ 

Remark 7.5. The above proof uses canonical form and therefore requires (l,cr) = (2,1). 
However, one can provide the alternative proof, similar to ones in subsection 2.3, which 
uses "precanonical" form (i.e. variables Qi,Q2)] now it is enough (I, a) = (I, or) >z (1,2). I 
leave details to the reader. 



7.3 Weakly-degenerate case 

Let us consider point x in which 

(7.21) I VV\ x C > C'\ log hr', (I, a)h(l, o-'). 

Then the same is true in its 7-vicinity with 7 = h s with some small 5 > provided 
C = C'(S) is large enough if (I, a) = (1 , cr') ; otherwise no restriction on C is needed. 
Thus all the arguments of the previous subsection remain valid but now for time T = e/i -1 
the shift with respect to £3 will be x (T provided |V||V| x ( and the shift with respect 
to (£2,^2) wn ^ De x Z 1 " 1 ^ provided (V^V^I x ( and therefore logarithmic uncertainty 
principle is fulfilled if 



(7.22) e > CC -1 M|log/i|- 

38 ' Actually it can be longer for small £3 but there is no benefit from this observation. 
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We need to assume also that 7 > Ce which is always the case provided fi < h s ° 1 with 
5 > 5. 

Then we can take 



(7.23) 




for I = 1, 
for / > 1 



and the contribution of this particular element ip 1 {x)q e {hD^) to the remainder estimate 
will be 0(p7 3 /i _2 T 1 ~ 1 ) (we assume that q is compactly supported). We can take always 
g = Cq\ and we estimated contribution of ^ 7 (x)(l — q e {hD^)) to the remainder estimate 
as 0(7 3 /i~ 2 ). So the contribution of ip-y(x) to the remainder estimate is 0(7 3 /i~ 2 ) anyway. 

Now let us consider its contribution to the approximation error. In our assumptions it 
does not exceed 

(7.24) C7r 2 7 3 + C Q ^h~ z d(e) x C7r 2 7 3 + C-f 3 h~ 3 (fih\ logh\) l+ ^C l \ \ogh\~ a 
and to keep it below C-y 3 h~ 2 we need to assume that 

21 + 1 , 2i-l . , , , 21-2cr-l 

(7.25) C > fi~h—\ \ogh\ ; 
in particular, for (I, a) = (1,2) 

(7.26) C>^^|log^|~' 
which is possible only for ji < h~^ \ \ogh\^. 

8 Spectral asymptotics 

This section finishes the analysis of 3-dimensional Schrodinger operator under condition 
(0.13). First of all, we have decomposition similar to one we had in two-dimensional case. 
Using this decomposition we prove that spectral asymptotics based on Ft — > /i~ 1 r(x T (t)r(ii'i/>)) 

with T = e/i _1 still holds with the remainder estimate 0(h~ 2 ) for \i < \i\ = f Ch~\ and thus 
prove theorem 0.3(h). 

Then we analyze the case \i > \i\ in two different settings: under non- degeneracy 
condition (0.11) we prove theorem 0.4(h), (hi) and without it we prove theorem 0.6. 
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8.1 Decomposition 

First of all, as I mentioned, we have decomposition, similar to (4.1)-(4.2): we reduce our 
operator to the family of 1-dimensional Schrodinger operators with respect to x 3 which are 
also /i _1 /i-pdos with respect to x 2 

(8.1) A n = hD 2 3 +r 2 n + W{x 2 ,x 3 ,fi~ 1 hD 2 ,fi~ 1 r n ), r n = ((2n + l)(ih)*, 
with x' = (xi,x 2 ) and arrive to decomposition 

(8.2) u(x,y,t) = T( ^u n (x',y',y 3 ,t)T n (x 1 )T n (y 1 ) JT" 1 
where u n (x', y', t) are Schwartz kernels of the propagators for A n . 



8.2 Proof of Theorem 0.3(ii) 

After above reduction , essentially we have 1-dimensional Schrodinger operators with respect 
to x 3 but with potentials which are /i _1 /i-pdos with respect to x 2 . We need to consider case 



(8.3) fl x = C7T3|log/i|-3 < ft < n^dhT* 

under assumption (/, a) = (1, 2); however, we consider general a > 2. Further, due to anal- 
ysis of sections 5,7 we need to consider only contribution to asymptotics and the remainder 
estimate of intermidiate zone 

(8.4) z' ^ { Ql < |£ 3 | < qi} with ei = c^- 1 

because we already know that the contribution of outer zone Z\ = {|£ 3 | > Qi} to the 
remainder is 0(h~ 2 ) and the contribution of the core zone 

(8.5) V={N<£i} 

to the remainder is 0{gifih~ 2 ) = 0{h~ 2 ). In definition (8.4) C\ is an arbitrarily large 
constant. 

Making 7-admissible partition with respect to x with 7 = h s we should consider two 

cases 

(8.6) ! |W|xC onsupp^', Ce[C,l], 
(8.6)2 \W\ < C = C\\ogh\- a on supp^' 
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where ip' is a 7-admissible element of the partition with respect to x. Further, taking 
in account the results of subsection 7.3, we can restrict ( from above by the right-hand 
expression in (7.26) (or (7.25)): 

(8.7) C < H*h*\ \ogh\-i; 

In what follows £ = ( in the case (8.6)2. 

Let us reduce operator to the canonical form and for any single n introduce temperate 
functions of (x 2 ,x 3 ,^ 2 ) 

(8.8) 7 = 7(n) = eC l \W + (2n + p = C(tC)*- 

We still keep e = C p~ l h\ log /i| and we consider now zone {7 > 7^} with j[ = C£ - s(/i| log h\)s 
(or equivalently {p > g[} with g[ = C((h\ \ogh\)s) chosen to have C _1 7 > e > C{p~ l h\ \ogh\)*. 

It is possible that g[ < Cp~ l ; in this case we redefine g[ as CpT x . 

First of all, using standard propagation of singularities one can see easily that 

(8.9) T f (Q'Q 3 u n ) = for T < \t\ < T x with T = Chp~ 2 \ log/i|, 

(rV+C) as (8.6)1, 



(8.10) Ti 



C _1 p as (8.6) : 



Here Q' = Q'(x2, x 3 , p, 1 hD 2 ) and Q 3 = q e {hD 3 ) are elements of the corresponding 7- 
admissible and ^-admissible partitions; q is supported in [— 1, — |] U [=, 1] and liere g = p. 
Note that 

(8.11) Tx > erV> ^0 < T x . 

Further, using standard propagation and rescaling technique one can see easily that 

(8.12) \F t ^ h -i T xAt)T'(Q f Q 3 u n )\ < CT P1 ^h~ 2 {^) s 

with T e [To, Ti] and To = p~ 2 h. Here and below x an d x are the same standard functions 
which we used before. 

Remark 8.1. (i) Due to ellipticity arguments, estimate 

(8.13) \F t ^ h -i TX T(t)r'(Q'Q 3 u n )\ < Ch s (l + g y s 

holds for q > C(p + |r|aV 

(ii) It also holds for p > C(g + |r|^) and g supported in {|£ 3 | < 1}. In this statement 
we do not need to assume that q vanishes in {|£ 3 | < |}. 

Therefore only elements with p x g provide non-negligible contribution. 
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Then inequality similar to (8.9) holds with Q' replaced by Q' which is an element of 
7-admissible partition of unity and with 7 replaced by 7; and now we can return to our 
original function u: 

(8.14) |F^ fc -i T (xT(t)r$Q 3 u))| < CTtffih- 2 ^)' x ^ 

where if) is an element of 7-admissible partition with respect to x and the factor g 2 /(ph) 
appears because there are no more than different numbers n which contribute. Here we do 
not have p anymore. 

It implies two important inequalities: 

(8.15) iF^h-wfartfQav))] < C Q fh~ 2 
and 

(8.16) \F t ^ T ((x Tl -Xf)(t)T(4,Q 3 u))\ < C(f?K-\^y 

as T = eyT 1 . 

The first of them implies that 

(8.17) |r(^g 3 e(.,.,o))-/i- 1 ^° (^ t ^r{xTM^Q^)^M< 

Cgfh^Ti 1 < Ctfh~ 2 . 

On the other hand, (8.16) implies that 

(8.18) {h-'j" ^(fe -Xf)(*)r(^g 3 «))d^| < C(ffh-\^) s . 

Then the same inequalities (8.17)-(8.18) hold with ip replaced by our target function ip, 
and with 7 replaced by 1. Further, summation with respect to partition of unity in £3 results 
in the inequalities with the same left-hand expressions and with right-hand expressions 
C(\\ogh\h~ 2 and Ch~ 2 (p 3 hy respectively because g[ > Cp^ 1 . Both these expressions are 
less than Ch~ 2 . 
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8.3 Proof of Theorem 0.3(ii). Conclusion 

Now we need to analyze zone {g[ > |^| > gi} if it is non-empty i.e. g[ > g± which is 
possible only in the case (8.6)i. In this zone we fix g = g[ and e = 7 = Ch\ \ogh\/ g± and 
consider just one element Q 3 . 

Further, one can study propagation of singularities with respect to x 2 , £2, £3 rather than 
with respect to x 3 . Then the same estimates (8.12),(8.14)-(8.18) hold with T\ = e( and 
To = Ch\ log/i|/(C^), Tq = Ch/((e) which are exactly the same as in subsection 8.2 with 
given fixed g and e. In the end of the day we get upper bounds Ch~ 2 again. 

Analysis in the zone {\C, 3 \ < £1} has been done in subsection 8.2. Summation of all of 
zones including {|£ 3 | > q~i} leads to the estimate 

(8.19) \T(t/)e(x,x,0)) - h" 1 f (F t ^ h -i T (xf(t)r^u)]dr\ < Ch~ 2 

J —00 

which in turn together with the standard result rescaled implies 

(8.20) I J tp(x) (e(x, x, 0) - £ w (x, 0)) dx\ < Ch' 2 

and one can see easily hat the approximation error does not exceed Ch~ 2 as well. 
► Theorem 0.3(h) is proven. 

Remark 8.2. The same arguments for \i > lead to the remainder estimate 0((/i/z)i/i~ 3 ) . 
This estimate is not the smoothness-related and cannot be improved without condition 
(0.11) even in the smooth case. 



8.4 Spectral asymptotics. h 3 < ^ < Ch 1 . Estimates. I 

As I mentioned, in this case we have essentially one- dimensional Schrodinger operator (or 
rather a family of them). 

We need to refine arguments of subsections 8.2, 8.3. Due to section 7 we need to consider 
zone Z\ = {|£ 3 | < g\ = C(fih\ log/i|)2 only. 

Let us introduce admissible functions 7 = 7( n ) and p = p( n ) of (X2, £3, £2) similar to 
those introduced in subsection 8.3: 

(8.21) 7 = eminjr/, \W + (2n + l)fih\ < ?%), \VW\ < p = ^( 7 )i 

Let us consider (7; ^)-admissible partition with respect to (^2, £3, £2; £3) and corresponding 
pdos Q'(x 2 ,x 3 , fi~ 1 hD 2 ) and Q 3 = q e {hD 3 ). Relations between p and g will be clarified 
later. 
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We assume that logarithmic uncertainty principle holds: 

(8.22) 07>C7i|log/i|, 7 2 > Cp' l h\\ogh\- 

this condition will be later weakened to fry > Ch,'-/ 2 > Cp~ 1 h. 

Note, that remark 8.1 (i) remains valid then and therefore we can assume that g < Cp; 
contribution of other partition elements will be negligible. 

Let us introduce in this zone (defined by (8.21) and (8.22)) mollification parameter 

(8.23) e = C r e _1 /i|log/i|. 
In this subsection we consider elements on which 

(8.24) \W + (2n+ l)/i/i| x p 2 = 0( 7 ); 

then remark 8. 1 (ii) remains valid as well and we can assume that g x p and we will not 
distinguish between these two parameters in this subsection. In particular, we need to 
consider only terms with numbers n' 

(8.25) \ n >- n \<c^-, g- P{n y 

Further, due to our standard propagation arguments T'(Q'Q 3 u n ) = for \t\ G [T ,Ti] 
with T = Chp~ 2 \ \ogh\ and T\ = e^p' 1 . 

Furthermore, due to standard propagation-rescaling arguments estimate (8.12) holds 
with To — Chp~ 2 ; here and below we use the same functions x? X as before. In turn this 
inequality implies inequality similar to (8.15): 

(8.26) \F^ h -^(xn'T{Q'Q3Un)) \ < C^W" 1 \r\ < eg 2 
In turn inequality (8.26) yields that 

(8.27) | Tr ((X^E(0)XT n , Q'Q 3 T n )- 

h- 1 J r^ ft -ir(xT(*)r(^Q , Q 3 «n)))dr| < 



where T is the reducing operator, E{r) is the spectral projector of the mollified and pulled 
back operator A, (., .) is an inner product in L 2 (M^), Tr' is an operator trace for operators in 
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L 2 (R 2 X ,), ip n v = ( < X~ 1 ij) < X(vT n ),T n ) for v = v(x 2 ,x 3 ) then ijj n = *jj n (x 2j x 3l p~ 1 hD 2 , p^hD^ 
and we will be more explicit later. 

Actually it was enough to prove (8.27) for g — 7 = 1, p > 1, h < 1, 

(8.28) |W+(2n+l)///i| < C, \VW\<C, 

\VW(z) - VW(z')\ < C|log|z-z'|| _1 on suppQ' 

and 

(8.29) \W+(2n + l)ph\ x 1 on suppQ'. 
Then, assuming instead that 

(8.30) |W + (2n + l)ph\ < Cq 2 , \VW\ < Cp 2 h, 

\VW(z) -VW(z')\ <Cp 2 ^ 1 -\\og\z- z'\y l on suppQ' 

and 

(8.31) \W+ (2n + l)ph\ x p 2 , onsuppQ' 

after rescaling x fc > Xk/l, h h/(gj), p 1— > /ry/g we get the previous case provided 

(8.32) p7 >Ch, p> p/7 

where for our choice of p, 7 the second condition is due to the first one and p > h~%. Also 
one can check easily that for our choice of p, 7 continuity condition is fulfilled after rescaling 
but one still need to check it in general. This analysis implies a little improvement of (8.27): 



Proposition 8.3. Estimate (8. 27) holds under assumptions (8.30) — (8. 

We also need a trivial 
Remark 8.4. Under conditions (8. 30), (8. 32) only (without (8.31)) both terms 
I Tr' ((%i/;E(0)%T n , Q'Q 3 T n )\ and [tr 1 J^F t ^ h -x T ^t)F(i; n Q'Q 3 u n )yjdT\ do not ex- 
ceed Cpj 3 ph~ 2 . In particular (8.27) holds as (pf x h. 
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8.5 Spectral asymptotics. h~s < fi < Ch~ l \. Calculations. I 

Further, ne can prove easily by the same methods as in section 5 of [Brlvr] the following 
Proposition 8.5. Under assumptions (8.30) — (8.32) estimate 

(8.33) \h~ l j (F t ^ h -i T (jc T (t)r(i; n Q'Q 3 u n )y\dT- 

V,/T 2 / 9 (-£! - W n ( X2 , x 3 , 6))^, ,3, 6, 60 ^ 2 d* 3 d& d6| < 

V| log T |- CT 

/to/ds cj n = ^n(^2, ^3,6, / i_1 6)Q'(a ; 2,a;3, 6)^3(6) ■ 

From this statement and proposition 8.3 we get immediately 
Corollary 8.6. Estimate 

(8.34) |Ty<(£^(0)yr n ,Q , Q 3 T n >- 

■^fih~ 2 J _ W»(z2,Z3,6))a>n(z2,a:3,6>£3) dx 2 dx 3 ^ 2 ^3| < 

holds. 

Remark 8.7. (i) With given £ = C7i| log(/i/p7)|^ _1 39 ** approximation error in (8.33) does 
not exceed C/i/i _1 7 2 | log7|" CT . 

(ii) Under conditions (8. 30), (8. 32) alone (8.34) holds for p7 x h. 

8.6 Spectral asymptotics. h~* < fi < Ch~ l . Intermediate zone. 
Estimates 

In this subsection we consider intermidiate zone 

(8.35) Z[ = { Ql = e (ph)? < |6| <Qi = Cmin((^|log/i|)3,l)} 

in which g 2 > e/ih and we need to take into account x q 2 j([ih) terms in the decomposition. 
This will compensate possible degeneracy. 

39 > Which is (8.23) rescaled. 
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We know that the contribution of (7; ^-admissible element of the partition to the re- 
mainder estimate (8.31) does not exceed Cp,h~ lr y 2 ; then 

(8.36) I WlV ((%i/;E(0)%T n ,Q'Q 3 T n )- 

n \ 

-^/i/i" 2 J 0^-£g - W n (x 2 , £3, ^ 2 ))^ n (x 2 , £3, 62, 6) dx 2 cb 3 d& d£ 3 



l-,2 w ^ _ r^-2 3 



because for no more than C^ 2 / (/x/i) numbers n we have p( n ) x g and in our case we should 
take -$(77) = 77 and p 2 x 7 because we need to satisfy condition (8.31); then for Q' replaced 
by 1-admissible pdo, the left-hand expression of (8.36) does not exceed Ch~ 2 . Further, 
after summation over ^-admissible partition with g = ||£ 3 | we get that for Q 3 supported in 
Z' the left-hand expression of (8.36) does not exceed 



(8.37) ChT 2 r g" 1 dg = Ch~ 2 \og{g x / Bl ) x Ch~ 2 C\ 

J 81 



^_idcf flog I log h\ ioip<h l \\ogh\ 1 , 

C 1 [ |log(/i/i)| for h-^hgh]- 1 < fi < Ch- 1 

which is a tiny bit larger than Ch~ 2 : it contains an extra factor £ _1 which does not exceed 
Clog I \ogh\. 

To get rid of this factor we need just to have |VVT| < (. Therefore let us introduce 
temperate function 7 by (8.21) with -#(77) = n\ log77|~ 2 . 

(i) Let us consider first elements with 7 < | \ogh\~ s and thus | VW\ < (; then p 2 = f C7I log7| 
C , ~f(. In turn, on such elements we can introduce temperate functions 

(8.38) £ = £ {n) = eC'lW + (2n + l)fih\, g = g {n) = < p 2 

with C = C an d we can make (£; ^-admissible (sub)partition in (x 2 , x 3 , £ 2 ; £3)- 

Then all the above arguments hold on these subelements and then for them the left-hand 
expression of (8.36) does not exceed C(£ 3 h~ 2 where g now is fixed and n are calculating 
accordingly. Therefore for original elements the left-hand expression of (8.36) expression 
does not exceed CC7 3 and the total contribution of such elements to the remainder estimate 
does not exceed (8.37) multiplied by £ and then we conclude that the contribution of all 
elements with 7 < | \ogh\~ s does not exceed Chr 2 . 



8 SPECTRAL ASYMPTOTICS 



71 



(ii) Now let us consider remaining elements. On each such element |VW| x ( = 
|log7|~ 2 G [CI]- Let us introduce again temperate functions £, g by (8.38) and consider 
corresponding (sub)partition in (x2, x$, £3)- 

On every subelement we can trace propagation until time T\ = eg which is larger than 
T = Chg~ 2 \ \ogh\ as long as g > C{h\ \ogh\)z which is always the case if 

(8.39) /I > /T^|log/i|§. 

But then for time T G [Ti, T*], T* = e^^ 1 the shift with respect x 3 does not exceed C(T; 
further, this shift is exactly of this magnitude provided |V||W| x £. Also, the shift with 
respect (0:2, £2) does not exceed C/i _1 CT; further, this shift is exactly of this magnitude 
provided |Vj_W| x (. 

Then logarithmic uncertainty principle holds as long as (g-e > Ch \ log/i| which holds for 
e = Ch\ log h\( ~ l g~ x and it does not exceed 7 under condition (8.39). Then the contribution 
of such subelement to the remainder does not exceed C g 2 t^^i~ x h~ 2 which in turn does not 
exceed C£ 3 \ log£>| -2 because obviously f? 2 7 _1 | < | log^|~ 2 . But then the total contribution 
of all such elements to the remainder does not exceed Ch~ 2 J £> -1 | log g\~ 2 dg < Ch~ 2 . 

On the other hand, contribution of every such element to approximation error does 
not exceed C h~ 3, y 3 gvarepsilon\ \ogh\~ 2 < Ch~ 2 -f 3 (~ 1 \ \ogh\~ 1 and their total contribution 
does not exceed Chr 2 \ \ogh\~\~ 2 . But this expression is less than Chr 2 . 

(iii) If condition (8.39) is violated one can apply the same arguments (with rescaling) 
but small parameter will be (h/ g 3 ) and terms, which were negligible before, now just do not 
exceed "correct estimate" multiplied by (h/g 3 ) s . However, this estimate of former negligible 
terms is sufficient to recover correct estimates. I leave all details to the reader. 

Finally, on each subelement the "correct Weyl expression" is 

(8.40) h- 1 f (F t „ h - w (xT(t)(r'Q'Q3u)))dT 

J —00 

with T = T = Chg~ 2 \ log/i| under condition (8.39) and with T = T = 
Chg~ 2 \ \og(h/ g 3 ) \ otherwise. One can notice easily that (8.40) is equivalent to 

(8.41) ^X^ -2 / ^("^3 ~ ^«( x 2,^3,6))^n(a;2,a;3,6,£3)^2^a;3^2^3 

n 

modulo Cl 3 \ \ogh\~ 2 and then the total contribution of these last error terms does not 
exceed Ch~ 2 . Therefore we have proven 

Proposition 8.8. The left-hand expression of (8.36) with Q' = I and Q3 supported in the 
intermediate zone (e.g. Q 3 = (x^ — Xgi)(hD 3 )) does not exceed Ch~ 2 . 
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8.7 Spectral asymptotics. Inner zone 

Now let us consider a true inner zone 

(8.42) Z" = {|£ 3 | <Qi = e min((^)3, l)}. 

Here we fix £> = £>i and will not refer to it as g. So far we do not distinguish between 
cases p < Co/i -1 and p > Coh~ l but in the latter case true inner zone covers classically 
permited zone completely. 

Let us introduce there functions 7( n ) and pr n ) by (8.21) but we take $(77) < 77 corre- 
sponding to declared smoothness which will be important now. We replace 7 by max(7, 70) 
with To defined by 

(8.43) 7o^(7o)^ =Ch 

to keep uncertainty principle (but not necessarily logarithmic uncertainty principle). 
Consider corresponding partitions of unity. 

(i) If on some element \W + (2n + l)ph\ x p 2 s then, according to corollary 8.6, con- 
tribution of each such element into remainder does not exceed Cph~ lr y 2 . Further, unless 
p < Qi, this contribution is negligible. 

(ii) Now consider elements with 

(8.44) \VW\ xC^pV 1 =0(7) "7" 1 - 

On each such element we define £, g functions of (x 2 , x 3 , £ 3 ) and n by (8.38) with £ 

— — 2 1 

defined by (8.40). We also replace £ by max(£, £q) where £$ = Chs(~3 is defined from 

(8.45) * .(C4)3 =C7i. 

On each subpartition element Let us consider corresponding subpartition we introduce 
mollifying parameter 

(8.46) e = Chg- 1 \\og(h/g£)\. 

Then for any subelement there exists only one number n with £>( n ) < g\. 

According to corrolary (8.6, contribution of each such subelement into remainder does 
not exceed Cph~ l £ 2 . Then taking the sum over all subelements, we derive thatthe contri- 
bution of the original element to the remainder estimate does not exceed 

/•min(ef/p 2 ,l) 

(8.47) CVi/rV / z- 1 dz = Cph~ l 1 2 \og{min{{g\/p 2 ,l) 1 /£ ). 

J loft 

This is a bit short of what we want: 
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Proposition 8.9. (i) Contribution of the partition element to the remainder does not exceed 

c^h~W- 

(ii) Contribution of the partition element to the difference between (8. 40) and (8. 41) 
does not exceed Cph~ lr -f 2 . 

(Hi) Contribution of the partition element to the mollification error does not exceed 
Cph^ 2 . 

Proof. Without any loss of the generality we can assume that p = 7 = 1. Really, we can 
always rescale x 1— > x/j, £3 1— > £3//?, £ 1— > £/j, g gj p, h 1— > h/(p r y), p 1— > p'j/p, e 1— > e/7, 

(i) To estimate the contribution of each subelement to the remainder by Cph~ x l 2 we 
applied proposition 8.3 with T\ = eg~ 1 £. However, we can increase its value a bit to 
T* = e^ _1 £-min(| log£|°", | \og(£h~z)\ a ); one can easily prove it by a long-range propagation 
result as in the proof of proposition (5.3). 

Then contribution of each partition subelement to the remainder does not exceed 

(8.48) Cph- l £ 2 (\\og£\- a + \\og(£h~i)\^. 

2 2 

This is true for £ > Chz and for £ x Chs as well because in the last case we can apply 
remark 8.7. 

(ii) Easy calculations show that contribution of each subelement to the difference be- 
tween (8.40) and (8.41) does not exceed Cph~ l £ 2 \ log£|" CT . 

Then summation of (8.48) with respect to subpartition brings estimate of the partition 
element to the remainder 

/.min( e ?,l) , v 

Cph~ x J 2 r 1 (j log^r CT + I log(^-3)|- CT j dx 2 dx 3 d& x ph~ l 

and the same estimate for the contribution of the partition element to the difference between 
(8.40) and (8.41). 

(iii) The contribution of each subelement to mollification error does not exceed Cph~ 2 g~ 1 £ 3 e\ log 
ph~ 1 £ 2 \ \ogh\~ 1 . Summation with respect to subpartition brings estimate of the contribu- 
tion of the element to the mollification error 

ph" 1 ] log /i| 1 J £~ x dx2 dx3 d^2 x ph~ l . 
Returning in (i)-(iii) to the original scale we get ph^ 1 1— > ph~ 1/ ~f 2 . □ 
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Since contribution of each partition element to the remainder does not exceed C\ih *7 2 , 
the whole contribution of Z" does not exceed 

(8.49) Cfih" 1 J \- 1 dx 2 dx 3 d£ 2 . 

This general answer translates into two different estimates: 

(i) Under nondegeneracy condition (0.11) we have 7x1 and the final estimate is Cuh~ x . 

(ii) In the general case the final estimate is C/i/i _1 7 _1 which is exactly second term in 
the right-hand expression of (0.27). So we had proven 

Proposition 8.10. (i) Under non- degeneracy condition contribution of inner zone Z" to 
the remainder estimate does not exceed C[ih~ x . 

(ii) In the general case for (I, a) -< (2, 0) contribution of inner zone Z" to the remainder 
estimate does not exceed 

(8.50) CyLhr 1 ^ 1 = C/i/i" 1 "^! log/i|-m 
with To defined by (8.43). 

Remark 8.11. One can easily prove statement (ii) for any (/, cr) ^ (/, a). Namely, assuming 
that 

(i) Coefficients have their m-th derivatives continuous with continuity modulus u(r]) = 
$(r])ri~ m with m = [(I, a)\ and that 

(ii) One of these derivatives is disjoint from (modified non-degeneracy condition) 
we, using induction by m and arguments of section 4.4 of [Ivrl] can recover estimate Cu~ l h. 

Then, using the same arguments we can drop the modified nondegeneracy condition 
and prove the same remainder estimate (8.50). 

8.8 Spectral asymptotics. Synthesis and clean-up 

Adding all three zones - outer (empty if u. > Coh~ 1 \ \ogh\~ 1 ), intermidiate (empty if u. > 
Coh~ l ) and inner, we conclude that the total remainder estimate is 0(h~ 2+tlh ) under 
non-degeneracy condition (0.11) and 0(hr 2 + /i/i" 1 " 1 ^ | log/il" 1 ^ in the general case. 

Now let us consider the total principal part. Note that the total contribution of Z\ = 
{l&l ^ f?i} to it is equal to 

( 8 - 51 ) Xj^" 1 / ^(e 3 2 + W / 'n(x 2 ,X3,6 + (2n + l)/i/i))x 

x(C3/0i)^n(x 2 , x 3 , f 2 , f 3 ) dx 2 dx 3 d£ 2 d£ 3 . 
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Easy calculations show that replacing in W n every copy of (2n + l)ph by W brings non- 
essential error which does not affect the accuracy of the final answer. Also one can see 
easily that replacing ip n by ip o <p we make an error lesser than h~ 2 + ph. 
So, adding the contribution of the exterior zone 

we get the final asnswer 

Y,^^ 1 I d(^ + V(x) + (2n+l) f ih)y i p(x)^)dx+ 

n>0 

E^' 1 / m/Qi)(o(£ + W(x) + (2n+l)»h)- 

n>0 ^ 

9(ti + V(x) + {2n + l)fjLhyji/>(x)y/g(xjdx 

and the first term is J S MW ip(x) dx while the second one is J £^^ijj(x) dx. 
► Theorems 0.4(h), (iii) and 0.3 are proven. 

9 Spectral asymptotics. Vanishing V or F. 

We want to apply rescaling method to get rid of conditions (0.13) V < — e and even (0.5) 
F>e . 

To get rid of condition (0.13) we introduce temperate functions 7(2;) = e|V(a;)| and 
p(x) = 7(2)2. Then rescaling ball -8(2,7(2;)) to 5(0,1) and dividing operator by p 2 we 
find ourselves again in our settings but with h, p replaced by h new = h/pj, p ncw = p^/p 
respectively. However we need a cut-off: in fact 7(2) = e| V(x)| + 70 with an appropriate 
parameter j . Modifications in some cases are needed. 

To get rid of condition (0.5) (while (0.13) holds) we introduce a temperate function 
7(2) = e|F(2)| + 7o with 7 = p~2. 

9.1 d = 2, vanishing V 

In this case we need to satisfy condition (0.13) only for theorem 0.1(i) and in calculations 
of correction term for p < \ogh\~ 1 and we do not need it for sure for p > h~ 1 \ logh\~ l . 
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Consider first the case of theorem O.l(ii) when fl 1 = h~s \ log | ~ a < ji < h" 1 . Since we 
do not need condition (0.13) for /i new ^new > 1, we will cut off 7 by picking up 70 = (fJ>h). 
Plugging /i new , 7 new to the remainder estimate (0.15) we see that the contribution of each 
element to the remainder does not exceed 

C/i-^-S + Cfji^hfy log(V7 f )|5- CT /r 2 7 3 "' 

and the total remaainder estimate will be then due to condition (0.11) 

C J At _1 /i7 _1 d 1 + Cin^tyh- 2 J |log(/i7-t)|^V" Z d 1 . 

The second integral is equal to the second term in the right-hand expression of (0.15) 
for I < 2; to handle the case / = 2 we need to notice that the second term in (0.15) 
is due to mollification error but in rescaled coordinates mollification parameter is e new = 
(/^new^newl log/inew|) 5 = (a 4 " 1 ^! l°g(^7~ 3 ) I ) 5 7 _1 an d in original coordinates it is then (/i _1 /i| log 
(fi~ l h\ log/i| 2 and thus mollification error does not exceed the second term in (0.15) anyway. 
Therefore 

(9.1) Theorem 0.1(ii) without condition (0.13) remains true with the remainder estimate 
C/x -1 /? -1 ! log(/i/i)| + R 2 where R 2 is the second term in the right-hand expression in (0.15). 



The analysis of the case of theorem 0.1(i) is even simpler: one need to pick 7 = /i -2 ; 
then on each element after rescaling either condition (0.13) is fulfilled or 7 x 70 and /i new x 1 
and condition (0.13) is not needed while h ncw x fi 3 h < 1; then one can prove easily that 

(9.2) Theorem 0.1(i) without condition (0.13) remains true with the remainder estimate 
C/j 1 ~ 1 h~ 1 \ log(fih)\ + Ri where R\ is the second term in the right-hand expression in (0.14). 

(9.3) Actually one can get rid of this logarithmic factor considering propagation of singu- 
larities which is going in direction orthogonal to VV with the speed x The actual 
velocity is f = fi" 1 ^ ■ (Fy l GVV where G = (g^ k ) and F = (F jk ) = (d Xh Vj - d Xj V k ); 
details in the smooth case can be found in subsection 6.5.4 of [Ivrl]. 

(9.4) Thus we can replace T\ = e/i by T\ = e//min(| log7|, | log(7/7o|) a in the Tauberian 
theorem provided (I, a) ^ (/, er). 

(9.5) Also, one can calculate "the second term" and to prove that it vanishes. 
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(9.6) After summation we get an extra term in the remainder estimate 

C^bT 1 J 7 - 1 (| log 7 |-< 7 + | log( 7 /7 |^' T ) dx x Cn^hT 1 . 

Therefore 

(9.7) Statements (i),(ti) of theorem (0.1) remain true with the same remainder estimates 
even without condition (0.13). 

9.2 d = 3, vanishing V 

For d = 3 we need condition without condition (0.13) only for fih < 1. Consider first 
Theorem 0.3 with \i < h~s. Then applying the standard scaling with p = 72 and 7o = /^~ 2 
we will find ourselves either in the case of condition (0.13) fulfilled or in the case of p ncw = 1 
when it is not needed. Then contribution of each partition element to the remainder is 
0(/i„ c 2 w ) = 0(/i -2 7 3 ) and then the total remainder is 0(h~ 2 ). Therefore 

(9.8) Theorem 0.3 holds without condition (0.13). Consider now h~s < \i < h^ 1 . Then 

(9.9) In non-degenerate case (under condition (0.11)) the remainder estimate was 0(h~ 2 ); 
plugging h ncw we get the contribution of such ball to the remainder estimate does not exceed 
Ch~ 2 p 2/ ~f 2 x Ch~ 2/ ~f 3 which results in the total contribution of all such balls 0(h~ 2 ). So, 

(9.10) Theorem 0.4(H) holds without condition (0.13). 

The same arguments show that in the total contribution of all such balls 0(h~ 2 ). So, 

(9.11) Theorem 0.4(i) holds without condition (0.13) provided I < 2. 

(9.12) In the general case we have an extra term R 3 = f C\ihM^h~ 2 \ log h]" 1 ^ in the re- 
mainder estimate. However, this term does not translates well: after we plug /x nGW , h new 
it produces Cph J + :i h~ 2 \ log(/i/p7)|~ I + 2T x 7 3_2 ('+ 2 ) and without condition (0.11) this greater 
than R 3 even if I — 1. 

However, we can introduce a different kind of scaling. Namely, introduce 7 defined by 

(9.13) 7 = mm{7] : \V\ < C$(r)), |W| < Cr]~H(r]). 



9 SPECTRAL ASYMPTOTICS. VANISHING V OR F. 



78 



I remind that (I, a) ^ (2,0); otherwise we can apply arguments of section 4.4 of [Ivrl] 
(see also Remark 8.11). In this case we can take p = ^\\og^\~^ and depending on 
which inequality transforms into equality we find ourselves after rescaling in two different 
situations: 

(i) If \V\ x 7 i |log7| _fJ we are in the degenerate situation with V ncw disjoint from 0. 
Plugging h new into Chr 2 we get Ch~ 2 p 2r y 2 which does not exceed Ch~ 2r y 3 . Plugging /i nC w, 
//new into Cph T +2h~ 2 \ log/i| _T +2 we get 



j n H a 2 l_ r, 

Cph'+^h log( — ) '+ 2 x p'+2 7 1+2 X7" 
PI 







(|log(-)| 


■ |log7l) 


V p7 



Cph'+ 2 h l [ | log(— )| • | log 7 1 ) 1+2 x 7 3 



which does not exceed C phJ+^h" 2 ] log/i| _r ?2 x 7 s . 

Summation over all these balls just eliminates factor 7 s . 

(ii) If I W| x 7 i_1 | log7|~°" then after rescaling we find ourselves in the non-degenerate 
situation with V new not necessarily disjoint from 0. However, we already examined this in 
(9.9)-(9.10) and found that contribution of the ball to the remainder estimate would not 
exceed Ch~ 2 w which in turn does not exceed Ch~ 2r y 3 and then the total contribution of such 
balls does not exceed Ch~ 2 . So, 



(9.14) Theorem 0.6 holds without condition (0.13). 



9.3 d = 2, vanishing F 

Assume that V is disjoint from but F vanishes and 
(9.15) \V(F 3 /V)\>e. 

Then we can introduce 7 = e\F\ + 70 and rescaling j(x)) to -6(0, 1) we find ourselves 
in frames of theorem 0.1 with /i new = p'j 2 , h new = h/ r y. Thus we need to take 7 = 
then for 7 x 70 we are in "normal magnetic field" situation (fi ncw x 1) when condition 
(0.5) is not needed and h new = pfih < 1 as long as p < h~ 2 . 

We assume that I — 2, a — 1. Consider first the case 1 < p < h~ x . Then the contribution 
of B(x, j(x)) to the remainder is 0(/i~ 1 /i" 1 7~ 1 ) and the total contribution does not exceed 
Cp^h- 1 J-f~ 3 dx x Ch- 1 due to (9.15). 

Consider now the case h~ l < p < h~ 2 . Let us assume that V x 1 (so o n = 
in condition (0.10)). Then contribution of B(x,j(x)) to the principal part is x h~ 2 ^ 2 
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and to the remainder 0(/i— l/i~ 1 7~ 1 ) as long as 7 < 71 = C(jjih)^ 1 . On the other hand, 
contribution of B(x, 7(2;)) to the principal part is and to remainder estimate is negligible. 
Then the total principal is x /i _2 7i x fi~ l h~ 3 and the remainder estimate is again. 

Consideration of extra terms due to the lesser smoothness is left to the reader. 
9.4 d = 3, vanishing F 

We assume that g^ k , F k have an extra smoothness 1 in comparison which is required in 
theorems we refer to. This assumption allows us to avoid problems straightening magnetic 
field. 

Consider first case fi < far 3. Then /i ncw < /i ne l as well and the contribution of 
B(x,7(x)) to the remainder estimate is 0(h~ 2 w ) = 0{h~ 2 ^ 2 ) and the remainder estimate 
does not exceed 

r cT^y 

(9.16) Ch~ 2 ^(x)- 1 dx< I Ch~ 2 log fi 

{ Ch~ 2 

with 

(9.17) i |VF|>e , 

(9.17)i ^|A jfe (VF)| >e 

j,k 

where A jfc (VF) denote 2x2 minors of matrix (VF) = (d Xj F k ). 

Note that (9.17)i implies that after rescaling condition (0.11) is fulfilled. 

Consider now case < fj, < h" 1 . Under condition (9.17)i contribution of B(x,^{x)) 
to the remainder is 0(h~ 2 w ) = 0(h~ 2r y 2 ) and then the total remainder is 0{h~ 2 log fi) and it 
is 0{h~ 2 ) under condition (9.17) 2 . Analysis in the case when (9.17)i fails is rather standard 
and boring. 

Finally, for h^ 1 < \l < h~ 2 magnitude of the remainder and the principal part strongly 
depend on the behavior of mes{x, F(x) < 77} as r] — > +0. In particular, in the generic case 

(9.17)3 |det(VF)|>e 

the remainder is 0(fi~ 2 h~ A ) while the principal part is of magnitude fi~ 3 h~ 6 . 



in the general case 
under condition (9.17)i 
under condition (9.17) 2 
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